Strange Quark Matter and Compact Stars 



F. Weber 



Department of Physics 
San Diego State University 
t^-". 5500 Campanile Drive, San Diego 

California 92182, USA 



O 
O 

(N 

o+ 

(N 
(N 

> 
in 
in 

o 
o 

o 

Of 

6 

c3 



February 2, 2008 



Abstract 



Astrophysicists distinguish between three different types of compact stars. These are white 
dwarfs, neutron stars, and black holes. The former contain matter in one of the densest forms 
found in the Universe which, together with the unprecedented progress in observational astronomy, 
make such stars superb astrophysical laboratories for a broad range of most striking physical 
phenomena. These range from nuclear processes on the stellar surface to processes in electron 
degenerate matter at subnuclear densities to boson condensates and the existence of new states 
of baryonic matter-like color superconducting quark matter-at supernuclear densities. More than 
that, according to the strange matter hypothesis strange quark matter could be more stable than 
nuclear matter, in which case neutron stars should be largely composed of pure quark matter 
possibly enveloped in thin nuclear crusts. Another remarkable implication of the hypothesis is 
the possible existence of a new class of white dwarfs. This article aims at giving an overview of 
all these striking physical possibilities, with an emphasis on the astrophysical phenomenology of 
strange quark matter. Possible observational signatures associated with the theoretically proposed 
states of matter inside compact stars are discussed as well. They will provide most valuable 
information about the phase diagram of superdense nuclear matter at high baryon number density 
but low temperature, which is not accessible to relativistic heavy ion collision experiments. 



1 



Contents 



ll Introducti 



ion 



2 Confined Hadronic Matter 



2.1 Effective Nuclear Field Theories] 



2.2 Non-relativistic Treatment, 



Primer on Quark Matter 



3.1 Models for the Equation of State 

3.2 Color- SuverconductivitA 



3.3 The Strange Quark Matter Hypothesis] 

3.4 Searches for Strange Quark Mattel 

3.5 Unusual Seismoaravhic Events] . . . . 

3.6 AMS and ECCC 



4 Relativistic Stellar Modelsl 



4.2 Stellar Structure Equations of Non-rotatina Stars] 



4.1 Particles in Curved Svace-Time 



4.3 Rotating Star Models] 



AA Kevler Freauenci 



4.5 Moment of Inertia of Rotating Compact Stars] 



5 Strangeness in Compact Starsl 



5.1 



Neutron Star, 



5.1.1 Huverom 



5.1.2 



K Meson Condensate 
Strange Quarks] . . . . 



5.2 



5.1.3 

5.1.4 H-dibaruonA 
Strange Stars] . . . 



5.2.1 Ceneral Provertie, 



sax J1SnS.J-Mf>A 



5.2.3 R.X J1856.5-375A 



5.2.4 The Neutron Star in 3C5t 



5.2.5 



5.2.6 



X-rau. Gamma-Rau Burst and SGR Associations] 
Rotational Instabilities] 



5.2.7 Surface Properties of Strange Stars] 



5.3 Strange Dwarfi 



6 Neutrino Emission and Stellar Cooling 

7 Signals of Quark Matter in Rotating Neutron Starsl 



7.1 Isolated Pulsar, 



7.2 Accreting Neutron Stars] 



8 Summary 



1 Introduction 



It is often stressed that there has never been a more exciting time in the overlapping areas of nuclear 
physics, particle physics and relativistic astrophysics than today. This comes at a time where new orbit- 
ing observatories such as the Hubble Space Telescope, Rossi X-ray Timing Explorer (RXTE), Chandra 
X-ray satellite, and the X-ray Multi Mirror Mission (XMM) have extended our vision tremendously, 
allowing us to see vistas with an unprecedented clarity and angular resolution that previously were only 
imagined, enabling astrophysicists for the first time ever to perform detailed studies of large samples of 
galactic and extragalactic objects. On the Earth, radio telescopes (e.g., Arecibo, Green Bank, Parkes, 
VLA) and instruments using adaptive optics and other revolutionary techniques have exceeded previous 
expectations of what can be accomplished from the ground. The gravitational wave detectors LIGO, 
LISA, VIRGO, and Geo-600 are opening up a window for the detection of gravitational waves emitted 
from compact stellar objects such as neutron stars and black holes. 

Neutron stars are dense, neutron-packed remnants of massive stars that blew apart in supernova 
explosions. They are typically about 20 kilometers across and spin rapidly, often making several hundred 
rotations per second. Many neutron stars form radio pulsars, emitting radio waves that appear from 
the Earth to pulse on and off like a lighthouse beacon as the star rotates at very high speeds. Neutron 
stars in X-ray binaries accrete material from a companion star and flare to life with a burst of X-rays. 
Measurements of radio pulsars and neutron stars in X-ray binaries comprise most of the neutron star 
observations. Improved data on isolated neutron stars (e.g. RX J1856. 5-3754, PSR 0205+6449) are now 
becoming available, and future investigations at gravitational wave observatories like LIGO and VIRGO 
will focus on neutron stars as major potential sources of gravitational waves. Depending on star mass 
and rotational frequency, gravity compresses the matter in the core regions of pulsars up to more than 
ten times the density of ordinary atomic nuclei, thus providing a high-pressure environment in which 
numerous subatomic particle processes compete with each other. The most spectacular ones stretch 
from the generation of hyperons and baryon resonances (£, A, 5, A) to quark (u, d, s) deconfinement to 
the formation of boson condensates (ir~, K~, H- matter) [H [21 El E] ■ There are theoretical suggestions 
of even more exotic processes inside neutron stars, such as the formation of absolutely stable strange 
quark matter j^J EJ Cj, a configuration of matter more stable than the most stable atomic nucleus, 
62 Ni. 1 In the latter event, neutron stars would be largely composed of strange quark matter [51 19*1 ITU] 
possibly enveloped in thin nuclear crusts [TTJ whose density is less than neutron drip. Another striking 
implication of the hypothesis is the possible existence of a new class of white dwarfs "1*2*1 H"""] . An 
overview of the conjectured composition of neutron stars is shown in Fig. [JJ Because of their complex 
interior structures, the very name neutron star is almost certainly a misnomer. Instead these objects 
should be coined nucleon stars, since relatively isospin symmetric nuclear matter-in equilibrium with 
condensed K~ mesons-may prevail in their interiors |14| . hyperon stars if hyperons (£, A, H, possibly 
in equilibrium with the A resonance) become populated in addition to the nucleons [13], quark-hybrid 
stars if the highly compressed matter in the centers of neutron stars should transform into u, d, s quark 
matter [T"""j, or strange stars if strange quark matter should be more stable than nuclear matter. The 
idea that quark matter may exist in the cores of neutron stars is not new but has already been suggested 
by several authors j"""""J "THJ, "TH1 1201 EH 1221 E3] For many years it has been thought that the deconfined 
phase of quarks and hadrons is strictly excluded from neutron stars. Theoretical studies, however, have 
shown that this was due to seemingly innocuous idealizations [2H 123 • Thus neutron stars may very well 
contain quark matter in their cores, which ought to be in a color superconducting state [2*6*1 [2*71 "2*51 [2*9*] . 
This fascinating possibility has renewed tremendous interest in the physics and astrophysics of quark 
matter. 

Of course, at present one does not know from experiment at what density the expected phase 

is common practice to compare the energy of strange quark matter to 56 Fe. The energy per particle of 56 Fe, 
however, comes in only third after 62 Ni and 58 Fe. 
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R ~ 1 km 



Figure 1: Competing structures and novel phases of subatomic matter predicted by theory 
to make their appearance in the cores (R £ 8 km) of neutron stars [2|. 



transition to quark matter occurs. Neither do lattice Quantum Chromodynamical simulations provide 
a conclusive guide yet. From simple geometrical considerations it follows that, for a characteristic 
nucleon radius of rjy ~ 1 fm, nuclei begin to touch each other at densities of ~ (47rr| i r/3)~ 1 ~ 0.24 fm~ 3 , 
which is less than twice the baryon number density of ordinary nuclear matter, p = 0.16 fm -3 (energy 
density eo = 140 MeV/fm 3 ). Depending on rotational frequency and stellar mass, such densities are 
easily surpassed in the cores of neutron stars so that gravity may have broken up the neutrons (n) 
and protons (p) in the centers of neutron stars into their constituents. Moreover, since the mass of the 
strange quark (s) is rather small, probably less than 100 MeV as indicated by the latest lattice results 
j30j . high-energetic up (it) and down (d) quarks may readily transform to strange quarks at about the 
same density at which unconfined up and down quarks appear. 

The phase diagram of quark matter, expected to be in a color superconducting phase, is very complex 
|2*d] \'27\. At asymptotic densities the ground state of QCD with a vanishing strange quark mass is the 
color-flavor locked (CFL) phase. This phase is electrically charge neutral without any need for electrons 
for a significant range of chemical potentials and strange quark masses pTT] . (Technically, there are 
no electrons only at zero temperature. At finite temperature the electron population is exponentially 
(exp(— A/T)) suppressed, where A denotes the superconducting gap.) If the strange quark mass is heavy 
enough to be ignored, then up and down quarks may pair in the two-flavor superconducting (2SC) phase. 
Other possible condensation patterns include the CFL-.fr phase |32] and the color-spin locked (CSL) 
phase [33] • The magnitude of the gap energy lies between ~ 50 and 100 MeV. Color superconductivity, 
which modifies the equation of state at the order (A//z) 2 level |3"4*ll35| . thus changes the volume energy 
by just a few percent. Such a small effect can be safely neglected in present determinations of models for 
the equation of state of neutron star matter and strange star matter. This is different for phenomena 
involving the cooling by neutrino emission, the pattern of the arrival times of supernova neutrinos, 
the evolution of neutron star magnetic fields, rotational (r-mode) instabilities, and glitches in rotation 
frequencies of pulsars (see Refs. j2Hl 123 HUH EZ1 EEl HUH HOI and references therein). Aside from neutron 
star properties, an additional test of color superconductivity may be provided by upcoming cosmic ray 
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space experiments such as AMS |H] and ECCO |42j . As shown in Ref. finite lumps of color-flavor 
locked strange quark matter, which should be present in cosmic rays if strange matter is the ground 
state of the strong interaction, turn out to be significantly more stable than strangelets jJUEH] without 
color-flavor locking for wide ranges of parameters. In addition, strangelets made of CFL strange matter 
obey a charge-mass relation that differs significantly from the charge-mass relation of strangelets made 
of ordinary strange quark matter ^0J|331 • This difference may allow an experimental test of CFL locking 
in strange quark matter [4*3] . 

In this review I will describe the current status of our understanding of the phases of superdense 
nuclear matter inside compact stars, putting special emphasis on the role of strange quark matter in 
astrophysics. This is accompanied by a discussion of possible observable signatures of the competing 
states of superdense matter in the cores of compact stars. These signatures will provide most valuable 
information about the phase diagram of superdense nuclear matter at high baryon number density but 
low temperature, which is not accessible to relativistic heavy ion collision experiments. 

The article is organized as follows. Section El discusses the properties and representative models 
for the equation of state of confined hadronic matter. This is followed by a brief primer on quark 
matter presented in section El Relativistic stellar models are discussed in section |U The possible role of 
strange quarks for compact stars and astrophysical phenomena associated with such stars are reviewed 
in section Neutrino emission from compact stars and their cooling behavior are discussed in section 
El Finally, possible astrophysical signals of quark matter in compact stars are reviewed in section 
followed by general concluding remarks provided in section |H1 

2 Confined Hadronic Matter 

The equation of state of neutron star matter below neutron drip, which occurs at densities around 
4 x 10 11 g/cm 3 , and at densities above neutron drip but below the saturation density of nuclear matter 
is relatively well known. This is to a less extent the case for the equation of state in the vicinity 
of the saturation density of normal nuclear matter. Finally the physical properties of matter at still 
higher densities are highly uncertain and the models derived for the equation of state of such matter 
differ considerably with respect to the functional dependence of pressure on density. This has its origin 
in various sources, which concern the many-body technique used to determination the equation of 
state, the model of the nucleon-nucleon interaction, the alteration of hadrons properties by immersion 
in dense matter, the fundamental constituents of neutron star matter (including phase transitions 
to meson condensates and quark matter), and the theoretical possibility that strange quark matter 
may be the true ground state of the strong interaction rather than nuclear matter. In the subsection 
below we introduce a collection of relativistic models for the equation of state which account for these 
uncertainties. Non-relativistic models for the equation of state will be studied in sectional 

2.1 Effective Nuclear Field Theories 

Up until the mid 1970s nearly all dense nuclear matter studies were based on non-relativistic potential 
models to describe the nucleon-nucleon interaction. The relativistic, field-theoretical approach to nuclear 
matter was pioneered primarily by Walecka and collaborators |JU1 \£7\ . The generalization of Walecka's 
model lagrangian to superdense neutron star matter has the following form [TJ 12 H3 EE] 

£ = E - m B ) ij B + l (d^od^o - mla 2 ) - ^F^F^ + ^m>^ + i (Pir-d^ - ml tt-tt) (1) 

B 

4r Zi t-, 4t?TZr 777 jr 
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+g P B^Bl^r-p^ B + -^-^b^t-G^b) - ^rn N b N (g aN a) 3 - ^c N {g aN af + J2i> L (i$ - m L ) ip L 

where B denotes baryons (p, n, E,A,H), L stands for leptons (e~,/i~), and of) = j^uj^, (j) = 7 M c^, 
u^ u = i[ji J, ) j l/ '\/2. The baryons are described as Dirac particles which interact via the exchange of 
a, u, it and p mesons. The a and uo mesons are responsible for nuclear binding while the p meson is 
required to obtain the correct value for the empirical symmetry energy The cubic and quartic a terms 
in Eq. (PJ) are necessary (at the relativistic mean-field level) to obtain the empirical incompressibility 
of nuclear matter |49j . The field equations for the baryon fields follow from as follows 01 

/ B 



Arris 

+ (g pB Y r-p M + J^-a^T-G^B + ^7 m 7 5 (^t.tt)^ . (2) 

The meson fields in (J2J) are obtained as solution of the following field equations, 

(d^df. + m^a = g^Bi'Bi'B - m N b N g aN (g aN a) 2 - c N g aN (g aN a) 3 , (3) 

B 

f B 

r ZTYIb 



+ ml) tt = £ — d " few^s) , (5) 

B Wl n 

d^G^ + m 2 p p u = Y^idpB^BT^B - -^<9 M (tPbtv^iPb)) , (6) 

B aTTIb 

with the field tensors F^ v and G^ u defined as F pu = d^oj v — d v uj^ and G pu = d^p u — d u p p . For neutron 
star matter, Eqs. (J2J) through © are to be solved subject to the conditions of electric charge neutrality 
and chemical equilibrium. The condition of electric charge neutrality reads 

£ 4(2J B + l)^f - E T% - Pm ®(P M - m M ) = , (7) 



B 



where Jb and q^ denote spin and electric charge number of a baryon B, respectively. The last term in 
(0) accounts for the electric charge carried by condensed bosons. The only mesons that may plausibly 
condense in neutron star matter are the tt~ or, alternatively, the more favored K~ OH] • Equation (|7jl 
constrains the Fermi momenta of baryons and leptons, kp B and kp L respectively. Leptons in neutron 
star matter are treated as free relativistic particles, 

(«V*0„ - m L) = . (8) 

The baryon and lepton Fermi momenta are further constrained by the condition of chemical equilibrium. 
Since neutron star matter is characterized by the existence of two conserved charges, electric charge 
and baryon charge, the chemical potential of an arbitrary baryon, B, created in a neutron star can be 
expressed in terms of two independent chemical potentials. Choosing p n and p e as the independent 
chemical potentials, one has 

H B = q B ^ n -qiv\ (9) 

with qs the baryon number of particle B. Since qs = <?/ = § for quark flavors f — u,d, s, the chemical 
potentials of quarks follow from © as 

P S = \p n -qf p e . (10) 
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Relativistic Greens functions constitute an elegant and powerful technique which allows one to derive 
from P |) to (|B |) a set of three coupled equations which are numerically tractable [2]. The first one of 
these equations is the Dyson equation which determines the two-point baryon Green function S B in 
matter, 

S B = S*-S* Z B S B . (11) 

The second equation determines the effective baryon-baryon scattering amplitude in matter, T BB ' , 
which is given by 

rpBB' _ yBB' _ yBB' _|_ f yBB j^BB' rpB'B' 

Here V and V ex denote the direct and exchange term of a given one-boson-exchange potential, which 
serves as an input, and the quantity A BB oc S B S B describes the propagation of baryons B and B' 
in intermediate scattering states. A popular and physically most suggestive choice for A is the so- 
called Brueckner propagator E2]- This propagator describes the propagation of two baryons in 
intermediate scattering states in terms of the full single-particle energy-momentum relation and, in 
addition, guarantees that these particles obey the Pauli principle too. The system of equations is closed 
by the expression for the self energy of a baryon in matter, S B , given by 

Z B = tJ2 / (Tr (T BB ' S b ')-T bb ' S b ') . (13) 

B' 

The one-boson-exchange potentials, V BB ' , sum the contributions arising from different types of mesons 
that exchanged among the baryons, 

(12\V BB '\?>A) = J2 5 4 (l,3)5 4 (2,4)r M (l,3)r M (2,4)A M (l,2), (14) 

M=cr,uj,iT,p,... 

where T M denote meson-nucleon vertices, and A M free meson propagators. Relativistic models for the 
equation of state of neutron star matter are obtained by solving Eqs. (fTTj) through (|T3*j) self-consistently 
in combination with Eqs. (JJJ) and Q for electric charge neutrality and chemical equilibrium. This has 
been accomplished for several different approximation schemes. With increasing level of complexity, 
these are the relativistic Hartree (RH), relativistic Hartree-Fock (RHF), and the relativistic Brueckner- 
Hartree-Fock (RBHF) approximations [2 j. The first two approximations are obtained by keeping only 
the Born term in the T-matrix equation (JUJ, that is, by setting T = V — V cx (and readjusting the 
parameters on the boson-exchange potential). The mass operator (fT3*j) is then given 

S B (P) = -i 9*b 9aB> /tS? eiW ° HO)^) - <WA> - q)S B {qj) ± . . . (15) 



B> 



where the T-matrix is replaced, according to Eq. (|14|1. by free meson propagators. This simplifies 
the solution of the many-body equations considerably. Because of this approximation, however, the 
coupling constants of the theory need to be adjusted to fit the properties of infinite nuclear matter, 
hypernuclear data and neutron star properties (see, for instance, Refs. [3 EH EJ ESI)- The five nuclear 
matter properties are the binding energy E/A, effective nucleon mass m^/mN, incompressibility K, 
and the symmetry energy a s at nuclear matter saturation density po (= 0.16 fm~ 3 ), 

E/A = -16.0 MeV, m* N /m N = 0.79, K ~ 225 MeV, a s = 32.5 MeV. (16) 

Of the five, the value for the incompressibility of nuclear matter carries the biggest uncertainty j^j. 
Its value is currently believed to lie in the range between about 220 and 250 MeV. In contrast to RH 
and RBHF, the RBHF approximation makes use of one-boson-exchange potentials whose parameters 
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e (MeV/fm 3 ) 

Figure 2: Models for the equation of state of neutron star matter computed for different 
compositions and many-body techniques described in the text (p,n denotes protons and 
neutrons; H, K~, Q stand for hyperons, K~ condensate, and quarks, respectively.) 



are adjusted to the properties of the deuteron and relativistic nucleon-nucleon scattering data. This 
approximation is therefore referred to as a parameter free. In passing we mention that in recent years 
a new class of effective field theories was developed which treat the meson-nucleon couplings density 
dependent. These field theories provide a very good description of the properties of nuclear matter, 
atomic nuclei as well as neutron stars |H2lEHllSZlEHlEni- We conclude this section with a brief discussion 
of the total energy density of the system follows from the stress-energy density tensor, T^ u , as j2] 

dC 

e = (T 00 ) = j£_ d^ x -g 00 C. (17) 

The pressure, and thus the equation of state, are obtained from Eq. (fTTjl as 

P = p 2 d/dp(e/p) . (18) 
Finally, the energy per particle is given in terms of the energy density and baryon number density as 

e = (E/A + m N )p. (19) 

Figure 121 shows several model equations of state based on RH and RBHF assuming different particle 
compositions of neutron star matter |2j. Of particular interest later (section |7j) will be the equation of 
state in Fig. El accounting for quark deconfinement pQ which, for this model, sets in at 230 ,MeV/fm 3 , 
which is less than two times the energy density of nuclear matter, eo = 140 MeV/fm 3 . Pure quark matter 
exists for densities greater than 950 MeV/fm 3 , which is around seven times eo- Of key importance for 
the possible occurrence of a quark-hadron phase in neutron stars is that pressure in the mixed phase 
of quarks and hadrons varies with density [21]. If this is not the case, hydrostatic equilibrium would 
strictly excluded the mixed phase from neutron stars. Table ^ shows the properties of nuclear matter 
computed for different many-body techniques as well as nuclear forces. In addition to the saturation 
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Table 1: Saturation properties of nuclear matter computed for numerous different nuclear 
forces and many-body techniques |2*1 I48| 15*4"]. 
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properties listed in (fTBj) this table also shows the symmetry energy density derivative, L, and the slope, 
y, of the saturation curve of isospin asymmetric nuclear matter, defined as L = Sp (da s /dp) po and 
y = — KdsfipoL)^, for several models, which vary considerably from one model to another. 

2.2 Non-relativistic Treatments 

The most frequently used non-relativistic treatments of dense nuclear matter studies are the hole-line 
expansion (Brueckner theory) jHOJ EI], the coupled cluster method, self-consistent Green functions, 
the variational approach (Monte Carlo techniques), the semiclassical Thomas- Fermi method [o^ I63j. 
and the density functional approach based on Skyrme effective interactions (for an overview of these 
methods and additional references, see Refs. (HI Ell)- Aside from the density functional approach, the 
starting point in each case is a phenomenological nucleon-nucleon interaction, which, in some cases, is 
supplemented with a three-nucleon interaction Vijk introduced to achieve the correct binding energy 
of nuclear matter at the empirical saturation density, po = 0.15 nucleons/fm 3 . Figure El shows two 
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Schroedinger-based models for the equation of state which are obtained for variational calculations 
based on the Urbana V\& two-nucleon interaction supplemented with the three-body interaction UVII 
(left curve) and TNI (right curve) three-nucleon interaction |3J ESI ESI- The Hamiltonian is thus of the 
form 

^ = E^v 2 + £^ + E Vijk- (20) 

i i<j i<j<k 

In the variational approach the Schroedinger equation 7i\^) = E\^f) is solved using a variational trial 
function |\P„), which is constructed from a symmetrized product of two-body correlation operators 
acting on an unperturbed ground state (Fermi gas wave function) |$), 

i*«) = (sn*«)i*>- ( 2i ) 

i<3 

The antisymmetrized Fermi-gas wave function is given by |$) = A(IIi<F|i)). The correlation operator 
contains variational parameters (14 for the UV14 and AV14 models jHH] and 18 for the more recent Aj.8 
model [HZj) which are varied to minimize the energy per baryon for a given density p, 

E(p) = min <*•'"'*■> , (22) 

which constitutes an upper bound to the ground-state energy of the system. Pressure and energy density 
of the stellar matter are obtained from relations (fTHj) and respectively. The new Thomas- Fermi 
approach of Myers and Swiatecki, TF96, is based on a Seyler-Blanchard potential generalized by the 
addition of one momentum dependent and one density dependent term [0*2"] . 

(23) 



Vn = Y(m) { kl =F Qa - i(l T C)(/j(f) 2 - 7^ + ')) ■ 



The upper (lower) sign in corresponds to nucleons with equal (unequal) isospin. The quantities 
kp , T (= kp o /2m), and po are the Fermi momentum, the Fermi energy and the saturation density of 
symmetric nuclear matter. The potential's radial dependence is described by a Yukawa-type interaction 
of the form 

1 p-ria/a 
Ana^ ri2/a 

Its strength depends both on the magnitude of the particles' relative momentum, p 12 , and on an average 
of the densities at the locations of the particles. The parameters £ and ( were introduced in order to 
achieve better agreement with asymmetric nuclear systems. The behavior of the optical potential is 
improved by the term a(2p/p ) 2 ^ 3 with the average density defined as p 2 / 3 = (p 2//3 + p^ 3 )/2, and pi 
and pi the densities of interacting neutron or protons at points 1 and 2. The seven free parameters of 
the theory are adjusted to the properties of finite nuclei, the parameters of the mass formula, and the 
behavior of the optical potential [U2] ■ The nuclear matter properties at saturation density obtained for 
TF96 are listed in table H 

3 Primer on Quark Matter 

3.1 Models for the Equation of State 

The field theory of quarks and gluons, Quantum Chromodynamics, is a non-Abelian gauge theory with 
SU(3) C as a gauge group. The QCD Lagrangian has the form 

C = r f {i lfl D» h - m f ) ^ - \F; v Fr , (25) 
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Table 2: Approximate masses, m/, and electric charge numbers, g^ 1 of quarks. 



Quark flavor (/) 


u 


d 


C 


s 


t 


6 


m/ (GeV) 
qf 


0.005 

T 3 


0.01 
i 

3 


1.5 

T 3 


0.1 

i 

3 


180 

+1 


4.7 
i 

3 



where ipj are the quark fields for each flavor / and m/ the current quark masses (see table |2J). In color 
space, the fields ipi are three-component columns with a = 1, 2, 3. The color gauge- covariant derivative 
]jv given by 

2& = tfo6^-«f [AiUG?, (26) 

where p s is the strong interaction coupling constant. The quantities Gf are the gluon fields with color 
indices i — 1, . . . , 8, and Aj the Gell-Mann SU(3) C matrices. The quantity F l is the gluon field tensor 
defined as 

*%> = d,Gl - d u Gl + g s fa Gj G k v , (27) 

where fyk are the SU(3) C structure constants. The equations of motion of the coupled quark and gluon 
fields derived from Eq. (|23|) are then as follows, 



(ij^-m f )r f = -gslA'i) (28) 




d»F; u + g s f ijk G^F* u = - g ,fy lv [^) ^. (29) 



Considerable efforts are made to solve the QCD equations of motion on the lattice. At present, however, 
such simulations do not provide a guide at finite baryon number density yet, and it is necessary to rely 
on non-perturbative QCD models for quark matter which incorporate the basic properties expected for 
QCD. Three different categories of models have emerged. These are (1) phenomenological models (MIT 
bag model) where quark masses are fixed and confinement is described in terms of a bag constant, and 
more advanced (2) dynamical models and (3) Dyson-Schwinger equation models where the properties 
of quark (matter) are determined self-consistently. The most widely used of these models is the MIT 
bag model [HBJ EH1 E3 ■ F° r this model, the pressure P % of the individual quarks and leptons contained 
in the bag is counterbalanced by the total external bag pressure P + B according to 

P + B = Y. pt > ( 30 ) 

/ 

while the total energy density of the quark flavors confined in the bag is given by 

e = J2t f + B. (31) 

/ 

The quantity B denotes the bag constant, and e are the contributions of the individual quark to the 
total energy density. The condition of electric charge neutrality among the quarks reads 

3£9?*S>-E *ft = °> ( 32 ) 

/ L 

where qf denotes the electric charge number of a quark of flavor /, listed in table El The contributions 
of each quark flavor to pressure, energy density and baryon number density are determined by the 
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10 2 fm 10 3 fm 10 4 fm 10 5 fm 10 6 fm 

Figure 3: Radii of quark bags [7T] . For masses less than 10 9 GeV the electrons (gray dots) 
are outside the quark bags (indicated by thick solid circles) and the core+electron system 
has a size of ~ 10 5 fm. For masses between 10 9 and 10 15 GeV the electrons are partially 
inside the core. For masses greater than 10 15 GeV all electrons are inside the core. 



thermodynamic potentials dVt f = —S f dT — P f dV — A^dp f , from which one obtains 

The quantity ^ denotes the chemical potential of quark flavor /, and rrif stands for its mass. The 
phase space factor Uf is equal to 2 (spin) x 3 (color) = 6. Chemical equilibrium among the quark flavors 
and the leptons is maintained by the following weak reactions, 

d <-» u + e~ + v e , s <-> u + e~ + v e , s c + e~ + v e , s + u <-> d + u , c + d <-> u + d . (34) 

Since neutron stars loose the neutrinos within the first few seconds after birth, the chemical potentials 
of neutrinos and anti-neutrinos obey [a u = fi u = and one obtains from the weak reactions (J34|) 

fjL d = H v + H e , fi c = fi u , fi = ft d = fi s . (35) 

The equation of state of relativistic quark matter at zero temperature made up of massless, non- 
interacting particles is readily calculated from (|33j). One obtains 

with Uf = 6. The equation of state of such matter is obtained from equations (j5Uj) and (JBTj) as 

P=(e-AB)/3. (37) 

From Eq. (|37jl one sees that the external pressure acting on a bag filled with quarks vanishes for e = 45. 
The mass contained inside the bag is given by M = f*edV = (4vr/3) eJ R 3 , which is the generic mass- 
radius relation of self-bound matter. The consequences of this relation for strange quark matter systems 
are illustrated in Fig. 01 The condition of electric charge neutrality of stellar quark matter, Eq. (|32|). 
leads to 

2p u - p d - p s = . (38) 
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Figure 4: Relative quark and lepton densities in quark-star matter as a function of mass 
density (from Ref. [72]). 



Since p u = p d = p s for massless quarks, one finds from Eq. ()30|) that for zero external pressure, P, 
the bag constant is related to the quark chemical potential as B = 3p 4 /ATT 2 . The energy per baryon 
number of quark matter follow as 



E 



B 



A p (p u + p d + p s )/3 p u 
with p the total baryon number density defined as 

/ 



4— = An 2 — 

p 3 



(39) 



(40) 



In the next step we turn to the determination of the equation of state if the charm and strange quarks 
are given their finite mass values listed in El In this case Eq. ()33)l leads for pressure, mass density, and 
baryon number density of quark matter to the following expressions, 



24vr 2 



3 1 + ,/l - zi 
2 J J 2 f 



z f 



-f 



z'f) ln 

2 f J 2 



z f 



(41) 
(42) 
(43) 



with Zi defined as Zf = rrif/pf. Ignoring the two most massive quark flavors, t and b, which are far too 
heavy to become populated in compact stars (see Fig. HJ), the condition of electric charge neutrality, 
expressed in Eq. fl3~2"|) . reads 

2{p u + p c ) - {p d + p s ) - 3(p e + ff) = . (44) 
Upon substituting (j4^|l into (|1H), this relation can be written as 

2(1 - {nVvf) (l + (1 - z 2 J) - (l + (1 - - (p e /p) 3 (l + (1 - = . (45) 
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Figure 5: Mass of quark stars at zero and 
finite temperature versus central star den- 
sity 
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Figure 6: Mass versus radius of the quark- 
star sequences shown in Fig. El [72] . 



An expression for the pressure of the system is obtained by substituting (J41j) into (|3(Jj) . This leads to 

,4 

\4 



P + B 



4vr 2 



(1 - (////z) 4 )^) + *(z a ) 



with $ defined as 



The total energy density follows from (J3"T|) as 



2 f 



+ ^zi In 



2 f 



(46) 



(47) 



(3P 



4B) + 4^5 E ^0 

/=s,c 



^ln- 



1 + 



1-4 



(48) 



The first term on the right-hand side of Eq. (|4Hj) represents the equation of state of a relativistic gas of 
massless quarks derived in Eq. ()37|) while the second term accounts for finite strange and charm quark 
masses. Finally, the total baryon number density of such quark matter is given by 



P= £s (C 1 - (^) 3 ) (l + (1 - ^ ) + (l + (1 - )) 



3tt 



(49) 



The relative quark-lepton composition of absolutely stable (I? 1 / 4 = 145 MeV) quark-star matter at zero 
temperature is shown in Fig. 0] All quark flavor states that become populated in such matter up to 
densities of 10 19 g/cm 3 are taken into account. Since the Coulomb interaction is so much stronger than 
gravity, quark star matter in the lowest energy state must be charge neutral to very high precision 9\. 
Therefore, any net positive quark charge must be balanced by a sufficiently large number of negatively 
charged quarks and leptons, which determines the lepton concentration shown in Fig. Because of 
their relatively large masses, the presence of charm quarks requires densities greater than 10 17 g/cm 3 in 
order to be present. Stellar sequences containing strange and charm quarks are shown in Figs. |3] and |3 
A stability analysis of these sequences against radial oscillations (see section I5~3)) . however, shows that 
only the strange star sequence is stable but not the charm star sequence. Finally, we mention that a 
value for the bag constant of B 1 / 4 = 145 MeV places the energy per baryon number of (non-interacting) 
strange quark matter at E/A = 829 MeV [7U], which corresponds to strange quark matter strongly 
bound with respect to 56 Fe whose energy per baryon is M( 56 Fe)c 2 /56 = 930.4 MeV, with M( 56 Fe) the 
mass of the 56 Fe atom. 
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Figure 7: Conjectured phase diagram for QCD [33]. For small rri^/A there is a direct 
transition from nuclear matter to CFL color superconducting quark matter. For large m? s / A 
there is an intermediate phase where condensation patterns like the CFL-K , CFL-fC + , 
CFL-7T , gCFL, 2SC, ISC, CSL, and LOFF phase (see text) may exist. Reprinted figure 
with permission from M. Alford, J. Phys. G 30 (2004) S441. Copyright 2004 by the Institute 
of Physics Publishing. 

3.2 Color- Superconductivity 

There has been much recent progress in our understanding of quark matter, culminating in the discovery 
that if quark matter exists it ought to be in a color superconducting state (2^1 E3 EHl EHj . This is made 
possible by the strong interaction among the quarks which is very attractive in some channels. Pairs 
of quarks are thus expected to form Cooper pairs very readily. Since pairs of quarks cannot be color- 
neutral, the resulting condensate will break the local color symmetry and form what is called a color 
superconductor. The phase diagram of such matter is expected to be very complex j2Hll2ZI, as can be 
seen from Figs. [7| and |H1 This is caused by the fact that quarks come in three different colors, different 
flavors, and different masses. Moreover, bulk matter is neutral with respect to both electric and color 
charge, and is in chemical equilibrium under the weak interaction processes that turn one quark flavor 
into another. To illustrate the condensation pattern briefly, we note the following pairing ansatz for the 
quark condensate |7ij . 



antisymmetric in Dirac indices, antisymmetric in color, and thus antisymmetric in flavor. The gap 
parameters A 1; A 2 and A 3 describe d-s, u-s and u-d quark Cooper pairs, respectively. The following 
pairing schemes have emerged. At asymptotic densities (m s — > or — > oo) the ground state of QCD 
with a vanishing strange quark mass is the color-flavor locked (CFL) phase (color-flavor locked quark 
pairing), in which all three quark flavors participate symmetrically. The gaps associated with this phase 
are 




A 3 ~ A 2 = Ai = A , 
and the quark condensates of the CFL phase are approximately of the form 



(51) 



(52) 
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Figure 8: Possible meson condensed phases in the neighborhood of the symmetric CFL 
state for a strange quark mass of m s = 150 MeV. fi is the quark chemical potential and 
fiQ the chemical potential for positive electric charge. Solid and dashed lines indicate first- 
and second-order transitions respectively [73]. Reprinted figure with permission from D. B. 
Kaplan and S. Reddy, Phys. Rev. D 65 (2002) 054042. Copyright 2002 by the American 
Physical Society. 



with color and flavor indices all running from 1 to 3. Since £ al3X £f a f b x = ^f a ^f b ~~ ^%^f a one sees that 
the condensate (p)2*|) involves Kronecker delta functions that link color and flavor indices. Hence the 
notion color-flavor locking. The CFL phase has been shown to be electrically neutral without any 
need for electrons for a significant range of chemical potentials and strange quark masses [3T]. If the 
strange quark mass is heavy enough to be ignored, then up and down quarks may pair in the two-flavor 
superconducting (2SC) phase. Other possible condensation patterns are CFL-.fr [22] > CFL-K + and 
CFL-tt '- [73], gCFL (gapless CFL phase) [73], ISC (single-flavor-pairing) [7J| [751 El , CSL (color-spin 
locked phase) [22], and the LOFF (crystalline pairing) [2TJ [Z7J [TH] phase, depending on m s , fi, and 
electric charge density. Calculations performed for massless up and down quarks and a very heavy 
strange quark mass (m s — ► oo) agree that the quarks prefer to pair in the two-flavor superconducting 
(2SC) phase where 

A 3 > , and A 2 = A 1 = . (53) 
In this case the pairing ansatz (|5U|) reduces to 

(^ 7 5<>ocAe a6 e^ 3 . (54) 

Here the resulting condensate picks a color direction (3 or blue in the example above), and creates 
a gap A at the Fermi surfaces of quarks with the other two out of three colors (red and green). The 
gapless CFL phase (gCFL) may prevail over the CFL and 2SC phases at intermediate values of m 2 s / 1 \i 
with gaps given obeying the relation (see Fig. EJ) 

A 3 >A 2 >Ai>0. (55) 

As shown in Fig. EH fc> r chemical potentials that are of astrophysical interest (see sections [S] through 
UJ, \i < 1000 MeV, the gap is between 50 and 100 MeV. The order of magnitude of this result agrees 
with calculations based on phenomenological effective interactions [23 UH] as well as with perturbative 
calculations for \i > 10 GeV [HO]. We also note that superconductivity modifies the equation of state at 



16 




50 100 150 
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Figure 9: Gap parameters Ai, A 2 , and A 3 as a function of m 2 / p for p = 500 MeV, in an 
NJL model p3j. There is a second order phase transition between the CFL phase and the 
gapless CFL (gCFL) phase at m 2 / p = 2 A. Reprinted figure with permission from M. Alford, 
C. Kouvaris, and K. Rajagopal, Phys. Rev. Lett. 92 (2004) 222001. Copyright 2004 by the 
American Physical Society. 

the order of (A/p) 2 [SUES!, which is even for such large gaps only a few percent of the bulk energy. Such 
small effects may be safely neglected in present determinations of models for the equation of state of 
quark-hybrid stars. There has been much recent work on how color superconductivity in neutron stars 
could affect their properties. (See Refs. [2H1 EE1 EH EE1 EH) and references therein.) These studies 
reveal that possible signatures include the cooling by neutrino emission, the pattern of the arrival times 
of supernova neutrinos, the evolution of neutron star magnetic fields, rotational stellar instabilities, and 
glitches in rotation frequencies. Several of these issues will be discussed in sections El to [7] 

3.3 The Strange Quark Matter Hypothesis 

The theoretical possibility that strange quark matter may constitute the true ground state of the strong 
interaction rather than 56 Fe was proposed by Bodmer [5], Witten jB], and Terazawa [7J. A schematic 
illustration of this so-called strange matter hypothesis is given in Fig. ^2 which compares the energy 
per baryon of 56 Fe and infinite nuclear matter with the energy per baryon of two and three-flavor quark 
matter. Three-flavor quark matter is always lower in energy than two-flavor quark matter due to the 
extra Fermi-well accessible to the strange quarks. Theoretical arguments indicate that the energy of 
three-flavor quark matter may be even smaller than 930 MeV in which case strange matter would be 
more stable than nuclear matter and atomic nuclei. This peculiar feature makes the strange matter 
hypothesis to one of the most startling speculations of modern physics, which, if correct, would have 
implications of fundamental importance for our understanding of the early universe, its evolution in 
time to the present day, astrophysical compact objects and laboratory physics, as summarized in table 
El [7(3 El IH3 Ell- I n the following we describe the possible absolute stability of strange quark matter 
for a gas of massless u,d,s quarks inside a confining bag at zero temperature [70]. For a massless 
quark-flavor /, the Fermi momentum, pp , equals the chemical potential, pj . The number densities, 
energy densities and corresponding pressures, therefore, follow from Eq. (}3l)|) as 

pf = (pff/ix 2 , e / = 3(/i / ) 4 /4 7 r 2 , P f = (m / ) 4 /4tt 2 = e / /3 . (56) 

For a gas of massless u and d quarks the condition of electric charge neutrality 2p u — p d = 0, which 
follows from Eq. (|52*|) as, requires that p d = 2p u . Hence the chemical potential of two-flavor quark 
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Figure 10: Superconducting gap A versus quark chemical potential p [ZB]- Reprinted figure 
with permission from D. B. Kaplan and S. Reddy, Phys. Rev. D 65 (2002) 054042. Copyright 
2002 by the American Physical Society. 



matter is given by p 2 = p u = p d /2 1 ^ 3 . The corresponding two- flavor quark pressure then follows as 
P 2 = p u -\-P d = (l + 2 4//3 )/i2/47r 2 = B. From the expressions for the total energy density, t 2 = 3P 2 + B = 
4B, and baryon number density, p 2 = (p u + p d )/3 = /i^/vr 2 , one then obtains for the energy per baryon 
of two-flavor quark matter, 



E 
A 



— = — = 934 MeV x B\{t , (57) 

p2 92 



where B^ = B^/Uh MeV. Values for B 1 ^ smaller than 145 MeV are excluded. Otherwise two- 
flavor quark matter would have a lower energy than 56 Fe, which then would be made up of up and down 
quarks, in contradiction to what is observed. In the next step we add massless strange quarks to the 
system. The three-flavor quark gas is electrically neutral for p u = p d = p s , i.e. p% = p u = p d = p s . For 
a fixed bag constant, the three-flavor quark gas should exert the same pressure as the two-flavor gas, 
that is P3 = P 2 . This implies for the chemical potentials ps = ((1 + 2 4 / 3 )/3) 1 / 4 /i2- Hence, the total 
baryon number in this case can be written as p 3 = pI/tt 2 = ((1 + 2 4 / 3 )/3) 3 / 4 p 2 - The energy per baryon 
is therefore given by 

~ A ^ = -^ = 829 MeVx<, (58) 
A 3 p 3 1.127p 2 

since e 3 = 3P% + B — 4B = e 2 . It thus follows that the energy per baryon of a massless non-interacting 
three-flavor quark gas is of order 100 MeV per baryon lower than for two-flavor quark matter. The 
difference arises from the fact that baryon number can be packed denser in three-flavor quark matter, 
P2/P3 = (3/(1 + 2 4//3 )) 3 / 4 ~ 0.89, due to the extra Fermi-well accessible to the strange quarks. The 
energy per baryon in a free gas of neutrons is equal to the neutron mass, E/A = 939.6 MeV. For an 56 Fe 
nucleus the energy per baryon is E/A = (56rriN — 56 x 8.8 MeV)/56 = 930 MeV, where = 938.9 
denotes the nucleon mass. Stability of two-flavor quark matter relative to neutrons thus corresponds 
to (E/A) 2 < m n , or B 1 ^ < 145.9 MeV (5 1 / 4 < 144.4 MeV for stability relative to 56 Fe). The stability 
argument is often turned around because one observes neutrons and 56 Fe in nature rather than two- 
flavor quark matter. Hence the bag constant must be larger than about 145 MeV. Bulk three-flavor 
quark matter is absolutely stable relative to 56 Fe for B 1 ^ < 162.8 MeV, metastable relative to a neutron 
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Figure 11: Comparison of the energy per baryon of 56 Fe and nuclear matter with the en- 
ergy per baryon of 2-flavor (u, d quarks) and 3-flavor (u, d, s quarks) strange quark matter. 
Theoretically the energy per baryon of strange quark matter may be below 930 MeV, which 
would render such matter more stable than nuclear matter. 



gas for B 1//4 < 164.4 MeV, and metastable relative to a gas of A particles for 5 1/4 < 195.2 MeV. These 
numbers are upper limits. A finite strange-quark mass as well as a non-zero strong coupling constant 
decrease the limits on B 1 ^ 4 |70[ ITU] . The presence of ordinary nuclei in nature is not in contradiction to 
the possible absolute stability of strange matter. The reason being that conversion of an atomic nucleus 
of baryon number A into a lump of strange quark matter requires the simultaneous transformation of 
roughly A up and down quarks into strange quarks. The probability for this to happen involves a weak 
transition oc G 2 ^ which makes nuclei with A £ 6 stable for more than 10 60 years. The conversion of 
very light nuclei into strange matter is determined by finite-size and shell effects which dominate over 
the volume energy of strange matter at small A values. An example for the mass formula of strange 
matter is [TU] H3| 143] 

— ~ (829 MeV + 351 MeV A' 2 ' 3 ) B] 1 ^ , (59) 

in which case strange matter becomes absolutely stable for A > 6. If quark matter is in the CFL phase, 
metastability or even absolute stability of strange quark matter may become more likely than hitherto 
thought since the binding energy from pairing of the quarks should reduce the energy of the system by a 
contribution proportional to A 2 Figure fT2l shows the energy per baryon for ordinary quark matter 
and CFL quark matter. For high A values a bulk value is approached, but for low A the finite-size 
contributions (surface tension and curvature) increase the energy per baryon significantly. The pairing 
contribution is on the order of 100 MeV per baryon for A »s 100 MeV for fixed values of the strange 
quark mass and bag constant. Another crucial difference between non-CFL and CFL quark matter 
is the equality of all quark Fermi momenta in CFL quark matter which leads to charge neutrality in 
bulk without any need for electrons [31] . This has most important consequences for the charge-to-mass 
ratios of strangelets. For non-CFL strangelets one has 

ZwO.lf — — V A for A < 10 3 , and Z & 8 ( — — V A 1 / 3 for A > 10 3 , (60) 

Vl50MeV/ V150MeV/ V ; 

while, in contrast to this, CFL strangelets have a charge-to-mass ratio of (13] 

Z « 0.3 ( ms ) A 2 / 3 . (61) 
V 150 MeV J v ; 
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Table 3: Strange matter phenomenology. 
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This difference may provide a test of color superconductivity by upcoming cosmic ray space experiments 
such as AMS jUJ and ECCO [12] (see section El and table HJ. 

3.4 Searches for Strange Quark Matter 

Experimental physicists have searched unsuccessfully for stable or quasistable strange matter systems 
over the past two decades. These searches fall in three main categories: (a) searches for strange matter 
(strange nuggets or strangelets) in cosmic rays, (b) searches for strange matter in samples of ordinary 
matter, and (c) attempts to produce strange matter at accelerators. An overview of these search 
experiments is given in table [3] The experiments searching for nuggets of strange matter which got 
stuck in terrestrial matter focus on objects whose masses can range from those of atomic nuclei to the 
upper limit of about 0.3 x IO -9 g. The latter carry a baryon number of A ~ 2 x 10 14 and have a radius 
of R ~ 10~ 8 cm. Strange nuggets heavier than 0.3 x 10~ 9 g will not be slowed down and stopped in the 
crust of the Earth. Finally nuggets of more than ~ 10 22 quarks (i.e. A ~ 10 21 ) would have too much 
momentum to be stopped by the encounter and thus would pass through the Earth. Such encounters 
could take the form of unusual meteorite events, Earth quakes, and peculiar particle tracks in ancient 
mica, in meteorites and in cosmic ray detectors [§8J. One distinguishing feature of unusual meteorite 
events caused by strange nuggets could be the apparent magnitude of —1.4 associated with a 20-gram 
nugget at a distance of 10 km, which would rival with that of the brightest star, Sirius. Another 
distinguishing feature could be the meteorite's velocity which is smaller than about 70 km/s for an 
ordinary meteorite bound to the solar system, but amounts to about 300 km/s for a strange meteorite 
of galactic origin. An upper limit on the flux of cosmic strange nuggets can be derived by assuming 
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Figure 12: Energy per baryon in MeV as a function of A for ordinary strangelets (dashed 
curves) and CFL strangelets (solid curves) for _B 1//4 in MeV as indicated, m s = 150 MeV, and 
A = 100 MeV jlHl- Reprinted figure with permission from J. Madsen, Phys. Rev. Lett. 87 
(2001) 172003. Copyright 2001 by the American Physical Society. 



that the galactic dark matter halo consists entirely of strange nuggets. The expected flux at the Earth 
is then on the order of 10 6 A' 1 V250 P24 cm~ 2 s" 1 sterad -1 , where P24 = p/(10~ 24 g/cm 3 ) and V250 is the 
speed in units of 250 km/s [10 . Experiments sensitive at this flux level or better have been able to rule 
out quark nuggets as being the dark matter for baryon numbers 10 8 < A < 10 25 |HTlll51j . This however 
does not rule out a low flux level either left over from the Big Bang or arising from collisions of strange 
stars. If the strange matter hypothesis is valid, one should indeed expect a significant background flux of 
nuggets from collisions of strange stars in binary systems, which are ultimately colliding because of the 
loss of angular momentum emitted from the binary system as gravitational radiation. If such collisions 
spread as little as 0.1 M® of strangelets with baryon number A (~ 10 3 ), a single collision will lead to 
a flux of 10 _6 A _1 t'25o cm" 2 s -1 sterad -1 [TU], assuming that the nuggets are spread homogeneously in 
a galactic halo of radius 10 kpc. This would lead to a concentration of nuggets in our galaxy of less 
than 10~ 8 nuggets/cm 3 , translating to a nugget concentration in terrestrial crust matter of at most 
10 9 cm -3 (mass density ~ 10~ 29 g/cm 3 ). Such a nugget density would corresponds to a concentration 
of nuggets per nucleon, Arrange /-Wnudeons, that is much less than 10~ 14 |152j . The upper limit on the 
concentration of strange nuggets per nucleon in terrestrial matter established experimentally by Briigger 
et al. |153j and Perillo Isaac et al. |141j is 10~ 14 , which falls short of the upper limits that follow from 
strange star-strange star collisions as well as he flux of strange nuggets in cosmic rays. Consequently 
the results of Briigger et al. |153j and Perillo Isaac et al. |141j shown in Fig. EH do not rule out the 
existence of strange matter. This figure shows the results for three samples, a heavy ion activation 
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Table 4: Past, present, and future search experiments for strange quark matter. 
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a AMS: Alpha Magnetic Spectrometer (scheduled for 2005-2008). 

b CRASH: Cosmic Ray And Strange Hadronic matter. 

c ECCO: Extremely- heavy Cosmic-ray Composition Observer. 

d 1MB: Irvine Michigan Brookhaven proton-decay detector (1980-1991). 

e JACEE: Japanese- American Cooperative Emulsion Chamber Experiment. 

f MACRO: Monopole, Astrophysics and Cosmic Ray Observatory (1989-2000). 

experiment mostly sensitive to light strangelets (A < 10 9 ) which, if present as cosmic rays, would 
have been absorbed into the Earth's atmosphere. This would be different for the Moon, which has no 
atmosphere. Since its surface has been exposed to cosmic rays for millions of years, the upper limit 
of the concentration of strange matter in the lunar soil can be used to deduce a limit for the flux of 
impinging strangelets. The lunar sample was collected from the top 0.5 to 1 cm surface, at the base of 
the Sculptured Hills, Station 8. The presence of high cosmic ray track densities in the sample suggests 
that the integrated lunar surface exposure age is about 100 Myr. Using the range of strange matter 
in normal matter suggested in jHH] , a limit on the flux of strangelets on the surface of the Moon was 
deduced in |141j and is shown in Fig. 

A limit on the total amount of strange matter in the Universe follows from the observed abundance 
of light isotopes. This is so because the strange nuggets formed in the Big Bang would have absorbed 
free neutrons which reduces the neutron to proton ratio, N n /N p . This effect in turn would lower the 
rate of production of the isotope 4 He, whose abundance is well known from observation. For a given 
mass of the strange nuggets, this constrains their total surface area. To be consistent with the missing 
dark matter, assumed to be strange quark matter, and the observed abundance of light isotopes, the 
primordial quark nuggets had to be made of more than ~ 10 23 quarks. According to what has been 
said above, quark nuggets that massive are not stopped by the Earth. 

As summarized in table EJ during the past years several experiments have been using high energy 
heavy-ion collisions to attempt to create strange quark matter (strangelets) in the laboratory. The 
detection of strangelets in relativistic heavy-ion collisions is conceptually rather simple, because of the 
strangelet's low Z/A ratio. Technically, however, such an attempt is very difficult because of several 
reasons. First of all one has to defeat the finite-number destabilizing effect of strange nuggets, that 
is, the number of free quarks produced in a nuclear collision must be sufficiently large enough such 
that surface and shell effects do not dominate over the volume energy of strange matter, destroying its 
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Figure 13: Experimental limits on the concentration of strangelets per nucleons, 
-^strange /-^Vnucieons, contained in three studied samples |141j . These are a meteorite, terrestrial 
nickel ore, and lunar soil. The results from Brugger et al. |153j obtained with an iron mete- 
orite are shown for comparison. Reprinted figure with permission from M. C. Perillo Isaac 
et al, Phys. Rev. Lett. 81 (1998) 2416. Copyright 1998 by the American Physical Society. 



possible absolute stability. Moreover, since the dense and hot matter is produced for ~ 10~ 22 s, there is 
no time to develop a net strangeness. However it is believed that strange nuggets will result from two 
types of simultaneous fluctuations that separate strange and antistrange quarks, and that also cool the 
nuggets so that they do not evaporate. Finally, we mention the huge multiplicity of particles produced in 
such collisions, which makes the particle identification rather cumbersome [154J. The NA52 (Newmass) 
experiment has been searching for long-lived strangelets as well as for antinuclei in Pb-Pb collisions at 
CERN SPS. No evidence for the production of long-lived charged strangelets has been observed. One 
very intriguing candidate for a strangelet of mass m = (7.4 ± 0.3) GeV, electric charge Z = —1, and 
laboratory lifetime r > 0.85 x 10~ 6 sec was detected in the data of the 1996 run |155| . This object, 
which could not be confirmed however, could have been made of Qu + 6d + 9s quarks, carrying a baryon 
number of A = 7, or 7u + 7d+ 10s quarks (A = 8). 

Evidence for the possible existence of strange matter in cosmic rays may come from Centauro cosmic 
ray events [H2J IHH1 E31 H23 [156 j. Such events are seen in mountain top emulsion chamber experiments. 
The typical energy of such an event is of order ~ 10 3 TeV, and the typical particle multiplicity is 
50 to 100 particles. Several intriguing Centauro events have been reported from a Brazilian- Japanese 
collaboration first several decades ago, where an interaction in the air 50 to about 500 meters above the 
detector gave rise to a large number of charged hadrons and zero or very few photons or electrons P-57J. 
In one particularly striking Centauro event, 49 hadrons were observed to interact in the detector but 
only one photon or electron. The typical transverse momentum (poorly determined) of about ~ 1 GeV 
seems to be larger than that for a typical event of the same energy. The striking feature is that there 
seem to be no photons produced in the primary interaction which makes the Centauro. This is unusual 
because in high-energy collisions, it mesons are always produced, and they decay into photons. A 
Centauro event is much like a nuclear fragmentation. If a nucleus were to fragment, then there would 
be many nucleons, and if the interaction which produced the fragmentation was sufficiently peripheral, 
there would be few pions. This possibility is ruled out because the typical transverse momentum is so 
large, and more important because a nucleus would not survive to such a great depth in the atmosphere. 
Being much more tightly bound together than an ordinary nucleus, a strangelet with a baryon number 
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Figure 14: Limits on the flux of strangelets impinging on the lunar surface |141j . Maximum 
cosmic flux refers to the cosmic flux of strangelets, assuming that all of the dark matter in 
the Universe is composed of strangelets. Reprinted figure with permission from M. C. Perillo 
Isaac et al, Phys. Rev. Lett. 81 (1998) 2416. Copyright 1998 by the American Physical 
Society. 



around A ~ 10 3 explains many of these unusual features. So it is conceivable that strangelets incident 
upon the top of the atmosphere or produced at the top of the atmosphere could survive to mountain 
altitude. It may have lost a significant amount of baryon number before getting to this depth however. 
A peripheral interaction might be sufficient to unbind it, since it certainly will not be so tightly bound 
with reduced baryon number. The problem with this explanation is that it does not explain the high 
transverse momenta. At transverse momenta of ~ 1 GeV one would expect final-state interactions 
to generate some pions, and therefore an electromagnetic component, which as mentioned above, is 
not observed |158j . Strange matter constitutes not the only explanation of Centauros. One alternative 
explanation would be that Centauro (and Ant i- Cent auro) events are manifestations of disoriented chiral 
condensates |159j . Another interpretation (of the Chacaltaya Centauro events) suggests that they are 
due to fragments of heavy primary cosmic rays. However the survival probability of heavy primary 
nuclei to this depth in the atmosphere appears much too low to account for the number of Centauros 
reported. 

Besides the peculiar Centauro events which may act as agents for strange matter, the high-energy 7- 
ray sources Cygnus X-3 and Hercules X-l may give additional independent support for the existence of 
strange matter. The reason is that air shower experiments on Her X-l and Cyg X-3 indicate that these 
air showers have a muon content typical of cosmic rays. This muon content is a surprising result. Typical 
cosmic rays at energies between 10 and 10 5 TeV are protons. To arrive from the direction of Cyg X-l 
or Her X-l, the particle must be electrically neutral. To survive the trip, it must be long-lived. The 
only known particle which can quantitatively satisfy this constraint is the photon. Photon air showers 
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however have only a small muon component. Photons shower by producing electron pairs. When only 
underground data was known, it was proposed that the most likely candidate for the initiating particle 
is a hadron, and in order for interstellar magnetic fields not to alter its direction relative to that of 
the source, the hadron-known in the literature as the cygnet-must be neutral. A natural candidate 
for the cygnet appears to be the H particle (see section 15.1.4)1 . the electrically neutral strangeness-2 
di-baryons with the quantum numbers of two lambdas [Z = 0, A = 2) proposed by Jaffe [86J. In the 
theory of hadrons composed of colored quarks and antiquarks, combinations other than the usual qqq 
and qq are allowed as long as they are color singlets. Jaffe found that a six-quark uuddss color-singlet 
state (H) might have strong enough color-magnetic binding to be stable against strong decay. That is, 
7tih could be less than the strong-decay threshold, twice the A (uds) mass, mn < 2m\o. Estimated 
lifetimes for H range from ~ 10~ 10 s for mn near the A°A° threshold to > 10 7 for light H particles near 
the nn threshold. The potentially long lifetimes raise the possibility that H particles may be present 
as components of existing neutral particle beams (e.g. E888 experiment listed in0J). To make the H 
long-lived enough, it is necessary to make the H have a mass below single weak decay stability. To 
generate a large enough flux of H particles, the source is assumed to be a strange star. Studies of 
Her X-l however seem to rule out this hypothesis, since studies of the correlation in arrival time with 
the known period of Her X-l give an upper limit of the particle mass of about 100 MeV. The source 
of radiation must be either due to anomalous interactions of photons or neutrinos, or from some exotic 
as yet undiscovered light-mass, almost stable particles. The problem with Cyg X-3 may be that it is 
accreting mass and thus has a crust, such that there is no exposed strange matter surface where small 
strangelets could be produced and subsequently accelerated electrodynamically to high energies into 
the atmosphere of the companion star where H particles were created via spallation reactions. 

Anomalously massive particles, which could be interpreted as strangelets, have been observed in a 
number of independent cosmic ray experiments jH5]. Two such anomalous events, which are consistent 
with electric charge values Z ~ 14 and baryon numbers A ~ 350, have been observed by a balloon-borne 
counter experiment devoted to the study of primary cosmic rays by Saito et al. |127j . A balloon-borne 
experiment carried out by the Italian/ Japanese CRASH (tabled} collaboration, however, could not 
confirm the existence of such objects in cosmic rays [128J. Evidence for the presence of strangelets in 
cosmic rays has also been pointed out by Shaulov |130t 1160] . This experiment, known as HADRON, 
was carried out at Tien-Shan Station between 1985 and 1993. Is is based on a combination of extensive 
air shower arrays and large emulsion chambers. The strangelet component in this experiment was 
estimated to be about 1 m~ 2 yr -1 . The data taken by HADRON indicate that some primary cosmic 
rays may contain non-nucleus components which generate extended air showers that contain both a 
large number of muons as well as very high energetic photons |160| . Another group of data, associated 
with the absorption of high energy photons in the atmosphere, suggests that cosmic rays may contain an 
unusual component with an absorption length a few times greater than for ordinary nuclei 160J. These 
features are nicely explained if one assumes that they are caused by stable or metastable strangelets 
160 . Besides that, the so-called Price event |161j with Z ~ 46 and A > 1000, regarded previously as 
a possible candidate for the magnetic monopole, turned out to be fully consistent with the Z/A ratio 
for strange quark matter [162 . Finally we mention an exotic track event with Z ~ 20 and A ~ 460 
observed in an emulsion chamber exposed to cosmic rays on a balloon has reported by Miyamura 129J. 
This exotic track event motivated the balloon-borne emulsion chamber experiment JACEE |132j and 
Concorde aircraft |163j experiments. JACEE was flown near the top of the atmosphere. At least two 
events have been observed which have been referred to as Ant i- Cent auros [133J. 

3.5 Unusual Seismographic Events 

As already described in section l3~4l De Riijula and Glashow speculated about the presence of lumps of 
stable strange matter, also referred to as strange nuggets or nuclearites, in the cosmic radiation |131j . 
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The seismic signals caused by these nuclearites passing through the Earth would be very different from 
the seismic signals caused by an earthquake or a nuclear explosion [SHI H31| . This follows from the rate 
of seismic energy produced by strange nuggets given by dE/dt = fapv 3 , where a is the nugget cross 
section, p the nominal Earth density, v the nugget speed, and / the fraction of nugget energy loss that 
results in seismic waves rather than other dissipation such as heat or breaking rock (5U|- Underground 
nuclear explosions have / ~ 0.01, chemical explosions / ~ 0.02. In contrast to this, strange nuggets 
with a mass of several tons (size of 10~ 3 cm) passing through the Earth would imply that / ~ 0.05. 
Anderson et al. [90J looked at more than a million records collected by the US Geological Survey between 
1990 and 1993 that were not associated with traditional seismic disturbances such as earthquakes. The 
seismic signature would be caused by the large ratio of the nuclearites speed, estimated to be around 
400 km/s. Strange nuggets might thus pass through the Earth at 40 times the speed of seismic waves. 
Most interestingly, Anderson et al. were able to single out two seismic events exhibiting this behavior. 
One event occurred on 22 October 1993, the other event occurred on 24 November 1993. In the first 
case, an unknown object seem to have entered the Earth off Antarctica and left it south of India. It 
was recorded at seven monitoring stations in India, Australia, Bolivia and Turkey. In the second case, 
an object seem to have entered in the South Pacific, south of Australia, and left the Earth 16.8 seconds 
later in the Ross Ice Shelf near the South Pole. This event was recorded at nine monitoring stations in 
Australia and Bolivia. The chord length between the entry and exit points of the 24 November 1993 
event is 4204 km so that the duration measured for this event, if caused by the passage of an object 
through the Earth, would imply a velocity for the hypothetical object of 250 km/s. The interpretation 
of the data as being caused by strange nuggets penetrating the Earth is backed by a Monte Carlo study 
that was used to identify the extent to which nuclearites could be detected by seismographic stations 
[9Til 1164] . The study showed that one would expect to detect as many as 25 4-ton nuclearite events per 
year if 4 ton strange nugget were to saturate the halo dark matter density. If 10% of the dark matter 
density were distributed in strange nuggets over the mass range from 0.25 to 100 tons one would expect 
about an event per year. Detection of a nuclearite would require at least six station sites to fix its 
impact time, location and velocity, and seismic detection of signals by at least seven stations is required 
in order to separate strange nuggets events from random spurious coincidences. 

3.6 AMS and ECCO 

As shown in Ref. > finite lumps of color-flavor locked strange quark matter, which should be present 
in cosmic rays if strange matter is the ground state of the strong interaction, turn out to be significantly 
more stable than strangelets without color-flavor locking for a wide range of parameters. In addition, 
strangelets made of CFL strange matter obey a charge- mass relation of Z/A oc A^ 1 ^ 3 , which differs 
significantly from the charge-mass relation of strangelets made of ordinary strange quark matter, as 
discussed in section I3~3l In the latter case, Z/A would be constant for small baryon numbers A, and 
Z/A oc A~ 2/3 for large A (TUJ E3 EH- This difference may allow an astrophysical test of CFL locking 
in strange quark matter The test in question is the upcoming cosmic ray experiment AMS-02 on 
the International Space Station scheduled jjJTJ. AMS-02 is a roughly a 1 m 2 sterad detector which will 
provide data from October 2005 for at least 3 years. AMS-02 will probe the dark matter content in 
various channels (anti-protons, anti-deuteron, e + , 7), cosmic rays and 7 astrophysics, and, as already 
mentioned, strangelets. The expected flux of strangelets of baryon number A < 6 x 10 6 at AMS-02 is 

m 

F ~ 5 x 10 5 (m 2 y sterad)" 1 x R_ A x M_ 2 x x t 7 , (62) 

where _R_4 is the number of strange star collisions in our galaxy per 10 4 years, M_2 is the mass of 
strangelets ejected per collision in units of 1O~ 2 M , Vjoo is the effective galactic volume in units of 
100 kpc 3 over which strangelets are distributed, and i? is the average confinement time in units of 
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10 7 years. All these factors are of order unity if strange matter is absolutely stable, though each with 
significant uncertainties |43j . 

Another intriguing instrument is ECCO (table HJ), whose primary goal will be the measurement of 
the abundances of the individual long-lived actinides (Th, U, Pu, Cm) in the galactic cosmic rays with 
excellent resolution and statistics. ECCO is a large array of passive glass-track-etch detectors to be 
exposed on orbit for at least three years (221 • The detectors will passively record the tracks of relativistic 
ultra-heavy galactic cosmic rays during exposure on orbit. After recovery, the detectors are calibrated, 
etched, and analyzed. ECCO is one of two instruments on the HNX (Heavy Nuclei eXplorer) spacecraft, 
which is under consideration as a Small Class Explorer Mission. The HNX mission is planned for launch 
in October 2005 into a 475-km circular orbit. Recovery is planned nominally for 3 years following launch. 



4 Relativistic Stellar Models 

4.1 Particles in Curved Space- Time 

Neutron stars are objects of highly compressed matter so that the geometry of spacetime is changed 
considerably from flat space. Thus models of such stars are to be constructed in the framework of 
Einstein's general theory of relativity combined with theories of superdense matter. The effects of 
curved space-time are included by coupling the energy-momentum density tensor for matter fields to 
Einstein's field equations. The generally covariant Lagrangian density is 

£ = C E + C G , (63) 

where the dynamics of particles are introduced through Ce and the gravitational Lagrangian density 
Cg is given by 

C G = g 1 / 2 R = g l ^g^R, l/ , (64) 
where g^ u and R^ v denote the metric tensor and the Ricci tensor, respectively. The latter is given by 

with the commas denoting derivatives with respect to space-time coordinates, e.g. , v = d/dx u etc. The 
Christoffel symbols T in (J65|) are defined as 

T % = \ 9 aX {9fi\, v + 9u\, » - 9^u, a) • (66) 
The connection between both branches of physics is provided by Einstein's field equations 

GT = BT - -g^R = 87cT^{e, P{e)) , (67) 

(/!, v = 0, 1, 2, 3) which couples the Einstein curvature tensor, G^, to the energy-momentum density 
tensor, T^ u , of the stellar matter. The quantities g^ v and R in (|67|) denote the metric tensor and the Ricci 
scalar (scalar curvature) _2j. The tensor T^ v contains the equation of state, P(e), on the stellar matter, 
discussed in sections |2] and El In general, Einstein's field equations and the many-body equations were 
to be solved simultaneously since the baryons and quarks move in curved spacetime whose geometry, 
determined by Einstein's field equations, is coupled to the total mass energy density of the matter. 
In the case of neutron stars, as for all astrophysical situations for which the long-range gravitational 
forces can be cleanly separated from the short-range forces, the deviation from flat spacetime over the 
length scale of the strong interaction, ~ 1 fm, is however practically zero up to the highest densities 
reached in the cores of such stars (some 10 15 g/cm 3 ). This is not to be confused with the global length 
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scale of neutron stars, ~ 10 km, for which M/R ~ 0.3, depending on the star's mass. That is to say, 
gravity curves spacetime only on a macroscopic scale but leaves it flat to a very good approximation 
on a microscopic scale. To achieve an appreciable curvature on a microscopic scale set by the strong 
interaction, mass densities greater than ~ 10 40 g/cm 3 would be necessary |165j ! This circumstance 
divides the construction of models of compact stars into two distinct problems. Firstly, the effects of 
the short-range nuclear forces on the properties of matter are described in a comoving proper reference 
frame (local inertial frame), where space-time is flat, by the parameters and laws of (special relativistic) 
many-body physics. Secondly, the coupling between the long-range gravitational field and the matter is 
then taken into account by solving Einstein's field equations for the gravitational field described by the 
general relativistic curvature of space-time, leading to the global structure of the stellar configuration. 



4.2 Stellar Structure Equations of Non-rotating Stars 

For many studies of neutron star properties it is sufficient to treat neutron star matter as a perfect 
fluid. The energy-momentum tensor of such a fluid is given by 

T*" = u^u u {e + P) + gVP, (68) 

where and u v are four-velocities defined as 

dx^ „ dx u , s 

They are the components of the macroscopic velocity of the stellar matter with respect to the actual 
coordinate system that is being used to derive the stellar equilibrium equations. The production of 
curvature by the star's mass is specified by Einstein's field equations, 

G^v = 8 7r T , where = -^g^R (70) 

is the Einstein tensor. The scalar curvature of spacetime R in Eq. (|70jh also known as Ricci scalar, 
follows from Eq. (J55J) as 

R = R^g flu . (71) 

Finally, we need to specify the metric of a non-rotating body in general relativity theory. Assuming 
spherical symmetry, the metric has the form 

ds 2 = _ e 2$(r) dt 2 + e 2A(r) ^2 + f 2 ^2 + r 2 ^ # 2 ? (?2) 

where $(r) and A(r) are radially varying metric functions. Introducing the covariant components of 
the metric tensor, 

g tt = - e 2 *« , g rr = e 2A « , g ee = r 2 , g H = r 2 sin 2 6 , (73) 
the non-vanishing Christoffel symbols of a spherically symmetric body are 

= e 2* W -2A W $ / (r) ; Vi = $ , (r) = A , (r) ^ r e & = r -l ) ^ = r -l ( T r g = _ re ~2 M r) 

-A cos 9 



sin 9 



r;, = - r sin 2 # e" 2 A « ,VL = - sin 9 cos 9 , (74) 



where primes denote differentiation with respect to the radial coordinate. From Eqs. (|68|h (|70|) and 
(174)1 one derives the structure equations of spherically symmetric neutron stars known as Tolman- 
Oppenheimer-Volkoff equations |Xf>f>| llf)7j 

dP _ e(r) m(r) (1 + P(r)/e(r)) (1 + 47rr 3 P(r)/m(r)) 
dr r 2 1 — 2m(r)/r 
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Figure 15: Neutron star mass versus radius for different equations of state. The broken 
horizontal lines refer to the masses of Vela X-l (1.86 ± O.16M ) [EE], Cyg X-2 (1.78 ± 
0.23 M Q ) [inn], and PSR 1913+16 (1.442 ± 0.003 M ) [IE]]. The line labeled 'TC denotes 
the average neutron star mass (1.350 ± 0.004 M ) derived by Thorsett and Chakrabarty 

m- 

Note that we use units for which the gravitational constant and velocity of light are G = c = 1 so that 
M & = 1.475 km. The boundary condition to ([75]) is P(r = 0) = P c = P(e c ), where e c denotes the energy 
density at the star's center, which constitutes an input parameter. The pressure is to be computed out 
to that radial distance where P(r = R) = which determines the star's radius R. The mass contained 
in a sphere of radius r (< R), denoted by m(r), follows as m(r) = 4ir J r dr' r' 2 e(r') . The star's total 
gravitational mass is thus given by M = m(R). 

Figure E3 shows the gravitational mass of non-rotating neutron stars as a function of stellar radius 
for several sample equations of state discussed in sections |2 and El Each star sequence is shown up to 
densities that are slightly larger than those of the maximum-mass star (indicated by tick marks) of each 
sequence. Stars beyond the mass peak are unstable against radial oscillations and thus cannot exist 
stably in nature. One sees that all equations of state are able to support neutron stars of canonical mass, 
M ~ 1.4 Mq. Neutron stars more massive than about 2 M on the other hand, are only supported by 
equations of state that exhibit a very stiff behavior at supernuclear densities and disfavor exotic (e.g., 
K~ mesons, quark matter) degrees of freedom. Knowledge of the maximum possible mass on neutron 
stars is of great importance for two reasons. Firstly, because the largest known neutron star mass 
imposes a lower bound on the maximum mass of a theoretical model. The current lower bound is about 
1.55 M & for neutron star Cyg X-2 [169., which does not exclude the existence of exotic phases of matter 
in the core of Cyg X-2. The situation could easily change if a future determination of the mass of this 
neutron star should result in a value that is close to its present upper limit of ~ 2M . The second 
reason is that the maximum mass is essential in order to identify black hole candidates |172t 1173] . For 
example, if the mass of a compact companion of an optical star is determined to exceed the maximum 
mass of a neutron star it must be a black hole. Since the maximum mass of stable neutron stars studied 
here is ~ 2.2 M , compact companions being more massive than that value would be black holes. 

4.3 Rotating Star Models 

The structure equations of rotating compact stars are considerably more complicated that those of non- 
rotating compact stars [2|. These complications have their cause in the rotational deformation, that 
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Figure 16: Features of a rotating compact star in general relativity. Indicated is the de- 
formation of the geometry of spacetime and the dragging of the local inertial frames. The 
latter rotate at a position dependent angular frequency uj(r,9,<p), which is to be calculated 
self-consistently from Einstein's field equations. The dragging frequencies inside three stellar 
configurations are shown in Fig. El 



is, a flattening at the pole accompanied with a radial blowup in the equatorial direction, which leads 
to a dependence of the star's metric on the polar coordinate, 9, in addition to the mere dependence 
on the radial coordinate, r. Secondly, rotation stabilizes a star against gravitational collapse. It can 
therefore carry more mass than it would be the case if the star would be non-rotating. Being more 
massive, however, means that the geometry of spacetime is changed too. This makes the metric functions 
associated with a rotating star depend on the star's rotational frequency. Finally, the general relativistic 
effect of the dragging of local inertial frames implies the occurrence of an additional non-diagonal term, 
g^, in the metric tensor g^ v '. This term imposes a self-consistency condition on the stellar structure 
equations, since the extent to which the local inertial frames are dragged along in the direction of 
the star's rotation, indicated schematically in Fig. is determined by the initially unknown stellar 
properties like mass and rotational frequency. The covariant components of the metric tensor of a 
rotating compact star are thus given by (21 1174] 

2 is i 2ib 2 2ib 2A 2a 2ib 

g tt = -e + e w uj , g t(f> = -e w u , g„ = e , gee = e * , = e w , (76) 

which leads for the line element to 

ds 2 = g^dx^dx v = - e 2u dt 2 + e 2 ^ (d<f) - to dt) 2 + e 2 » d6 2 + e 2X dr 2 . (77) 

Here each metric function, i.e. v, if), \x and A, as well as the angular velocities of the local inertial frames, 
u, depend on the radial coordinate r and polar angle 9, and implicitly on the star's angular velocity Q. 
Of particular interest is the relative angular frame dragging frequency, ou, defined as 

u)(r,9,Q) = -u(r,9,n) , (78) 

which is the angular velocity of the star, Q, relative to the angular velocity of a local inertial frame, u. 
It is this frequency that is of relevance when discussing the rotational flow of the fluid inside the star, 
since the magnitude of the centrifugal force acting on a fluid element is governed-in general relativity 
as well as in Newtonian gravitational theory-by the rate of rotation of the fluid element relative to a 
local inertial frame |175j . In contrast to Newtonian theory, however, the inertial frames inside (and 
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Figure 17: Dragging of the local inertial frames inside rotating neutron stars in equatorial 
direction. Pk and M denote Kepler period, defined in Eq. (jHSJ), and gravitational mass. 
(Fig. from Ref. 0) 

outside) a general relativistic fluid are not at rest with respect to the distant stars, as pointed out just 
above. Rather, the local inertial frames are dragged along by the rotating fluid. This effect can be quite 
strong, as shown in Fig. El For a heavy neutron star rotating at its Kepler frequency, one sees that 
uj/VL varies typically between about 15% at the surface and 60% at the center, where the mass density 
is highest. 

4.4 Kepler Frequency 

No simple stability criteria are known for rapidly rotating stellar configurations in general relativity. 
However an absolute upper limit on stable neutron star rotation is set by the Kepler frequency Qk, 
which is the maximum frequency a star can have before mass loss (mass shedding) at the equator sets 
in. In classical mechanics, the expression for the Kepler frequency, determined by the equality between 



the centrifugal force and gravity, is readily obtained as Qk = yM/R 3 . In order to derive its general 
relativistic counterpart, one applies the extremal principle to the circular orbit of a point mass rotating 
at the star's equator. Since r = 6 = const for a point mass there one has dr = d6 = 0. The line element 
(177)1 then reduces to 





(79) 



Substituting this expression into J = ds, where si and S2 refer to points located at that particular 
orbit for which J becomes extremal, gives 




(80) 




(81) 



(82) 



31 



0.3 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



M/M 



o 



Figure 18: Onset of mass shedding from rapidly spinning neutron stars, computed for a 
collection of equations of state (Fig. from Ref. [2]). The Kepler period is defined in Eq. (|83|). 



It constitutes a simple quadratic equation for the orbital velocity V of a particle at the star's equator. 
Solving (|81|) for Q = fl^ gives the fluid's general relativistic Kepler frequency in terms of V, the metric 
functions v and ip, and the frame dragging frequency u, each quantity being a complicated function of 
all the other quantities. In this manner is given by j2] 



which is to be evaluated at the star's equator (primes denote radial derivatives). Figure fTSI shows Pk 
as a function of rotating star mass. The rectangle indicates both the approximate range of observed 
neutron star masses as well as the observed rotational periods which, currently, are P > 1.6 ms. One 
sees that all pulsars so far observed rotate below the mass shedding frequency and so can be interpreted 
as rotating neutron stars. Half-millisecond periods or even smaller ones are excluded for neutron stars 
of mass 1.4 M & |2~l ll74|ITTo r ] . The situation appears to be very different for stars made up of self-bound 
strange quark matter, the so-called strange stars which will be introduced in section 15.21 Such stars 
can withstand stable rotation against mass shedding down to rotational periods in the half-millisecond 
regime or even below |152j . Consequently, the possible future discovery of a single sub-millisecond 
pulsar spinning at say 0.5 ms could give a strong hint that strange stars actually exist, and that the 
deconfined self-bound phase of three-flavor strange quark matter is in fact the true ground state of the 
strong interaction rather than nuclear matter. Strange stars of a canonical pulsar mass around 1.4 M & 
have Kepler periods in the range of 0.55 ms £ Pk ~ 0.8 ms, depending on the thickness of the nuclear 
curst and the bag constant [H1H2]- This ran g e is to be compared with Pk ~ 1 ms obtained for standard 
(i.e., no phase transition) neutron stars of the same mass. Phase transitions in neutron stars, however, 
may lower this value down to Kepler periods typical of strange stars |177j . 




(83) 
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4.5 Moment of Inertia of Rotating Compact Stars 

To derive the expression for the moment of inertia of a rotationally deformed, axisymmetric star in 
hydrostatic equilibrium, we start from the following expression, 

I (A, fi) = - JjLr d6 d<f) T/(r, 9, 0, SI) y/-g(r,e,(/),n) . (84) 

We assume stationary rotation, which is well justified for our investigations. The quantity A denotes 
an axially symmetric region in the interior of a stellar body where all matter is rotating with the same 
angular velocity Q, and y/—g = e x +^+ u +v } _ The component of the energy-momentum tensor is given 
by 

= (t + P)u <t> u t . (85) 

Let us focus next on the determination of the fluid's four- velocity, u K = (u f , u r , u 9 , v?). From the general 
normalization relation u K u K = — 1 one readily derives 

-1 = (u t ) 2 g tt + 2u t u' t> g^+(u' t> ) 2 g H . (86) 

This relation can be rewritten by noticing that 

= VLu 1 , (87) 

which extremizes the total mass energy of the stationary stellar fluid subject to the constraint that the 
angular momentum about the star's symmetry axis, J z , and its baryon number, A, remain fixed [178J. 
Substituting (JBTj) into (|8Bjl leads to 

-1 = (u t ) 2 (g tt + 2g t ^Q + g^Q 2 ), (88) 

which can be solved for it*. This gives 

u* = (-(g u + 2 gut, O + g H Q 2 )) . (89) 



Replacing g a , gt<f> and g^ with the expressions given in (|7H|) and rearranging terms leads for (|%9*j) to 

e~ v (l - (uj - ty 2 e^- 2 ^ 1 . (90) 



Last but not least we need an expression for of Eq. (|85|) in terms of the star's metric functions. To 
accomplish this we write as = g^ R u K = g ( f > tu t + g^vr . Upon substituting the expressions for g^ t 
and from ()76|) into this relation, we arrive at 

U(f> = (Q -u) e 2lp u l . (91) 

Substituting the four-velocities (J5U|l and (J9~T]) into (J85j) gives the required expression for the energy- 
momentum tensor, 

T/ = (t + P){n-uj)e 2lp (e 2u -{uj-nfe^Y 1 . (92) 



Finally, inserting this expression into (J84j) leads for the moment of inertia of a rotationally deformed 
star to 

(93) 



I (Q) = 2vr / d9 / dr 
Jo Jo 



Relativistic corrections to the Newtonian expression for J, which for a sphere of uniform mass density, 
e(r) = const, is given by I = |Mi? 2 , come from the dragging of local inertial frames ((Q — uj)/VL < 1) 
and the curvature of space. 
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Figure 19: Mass versus radius of neutron stars, white dwarfs, and planets. Stars located between the 
lightest neutron star (marked '2') and the heaviest white dwarf (marked '3') are unstable against radial 
oscillations and thus cannot exist stably in nature. If strange matter would be more stable than nuclear 
matter, an enormous region in the M-R plane void of compact objects could be populated with strange 
matter configurations (see Fig. l2*5Jl . 

5 Strangeness in Compact Stars 

5.1 Neutron Stars 

Physicists know of three types of compact stars. These are black holes, neutron stars, and white 
dwarfs. Neutron stars and white dwarfs are in hydrostatic equilibrium so that at each point inside 
the star gravity is balanced by the degenerate particle pressure, as described mathematically by the 
Tolman-Oppenheimer-Volkoff equation (|75j) . These stars, therefore, exhibit the generic mass-radius 
relationship shown in Fig. EH Depending on composition the maximum neutron star mass (marked 
'1') lies somewhere between ~ 1.5 and 2.5 M® (see also Fig. IT5jl . The minimum mass of neutron 
stars (marked '2') is around 0.1 M®. These stars are considerably less dense (~ O.leo) and thus much 
bigger (R £j 300 km) than the neutron stars of canonical mass, which is ~ 1.4 M . White dwarfs at 
the Chandrasekhar mass limit (marked '3') have densities around 10 9 g/cm 3 , which is five orders of 
magnitude smaller than the typical densities encountered in neutron stars. 

5.1.1 Hyperons 

Model calculations indicate that only in the most primitive conception, a neutron star is made of only 
neutrons. At a more accurate representation, a neutron star may contain neutrons (n) and protons (p) 
whose charge is balanced by electrons (e~) and muons (/i~), strangeness-carrying hyperons (E,A,S), 
meson condensates (K~ or p~~), u, d, s quarks, or possibly H-dibaryons. The particle composition is 
determined by the conditions of electric charge neutrality and chemical equilibrium as well as the in- 
medium properties of the constituents calculated for a given microscopic many-body theory (cf. sections 
|2]and|2I). In general, the population of negatively charged hadronic states is favored over the population 
of positively charged hadronic states since the negative charge carried by the hadrons can replace high- 
energy electrons, by means of which a lower energy state is reached. Aside from electric charge, the 
isospin orientation of the possible constituents is of key importance for the population too. The reason 
is that neutron star matter constitutes a highly excited state of matter relative to isospin symmetric 
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Figure 20: Effective kaon mass in nuclear |183j and neutron star |184j matter. (Fig. from Ref. [2].) 



nuclear matter. Hence, as soon as there are new hadronic degrees of freedom accessible to neutron star 
matter which allow such matter to become more isospin symmetric, it will make use of them. 



5.1.2 K Meson Condensate 

The condensation of K~ mesons in neutron stars is initiated by the reaction 

e~ — > K~ + v . 



(94) 



If this reaction becomes possible in a neutron star, it is energetically advantageous for the star to 
replace the fermionic electrons with the bosonic K~ mesons. Whether or not this happens depends on 
the behavior of the K~ mass in neutron star matter. Experiments which shed light on the properties 
of the K~ in nuclear matter have been performed with the Kaon Spectrometer (KaoS) and the FOPI 
detector at the heavy-ion synchrotron SIS at GSI |179t I180| I181| I182j . An analysis of the early K~ 
kinetic energy spectra extracted from Ni+Ni collisions [179 J showed that the attraction from nuclear 



matter would bring the K 
relation 



mass down to m 



200 MeV at p ~ 3 po- For neutron-rich matter, the 



K 



(p) ~ m K - 1 



0.2 — 



(95) 



was established |185| I186| I187j . with rriK = 495 MeV the K~ vacuum mass. Values around m 
200 MeV lie in the vicinity of the electron chemical potential, p e , in neutron star matter [21 so that 
the threshold condition for the onset of K~ condensation, p e = m* K , which follows from Eq. (j94j) . could 
be fulfilled in the centers of neutron stars. The situation is illustrated graphically in Fig. Equation 
(jlSJ) is followed by 

n + e~^p + K~ + u, (96) 

with the neutrinos leaving the star. By this conversion the nucleons in the cores of newly formed 
neutron stars can become half neutrons and half protons, which lowers the energy per baryon of the 
matter 188.. The relative isospin symmetric composition achieved in this way resembles the one of 
atomic nuclei, which are made up of roughly equal numbers of neutrons and protons. Neutron stars are 
therefore referred to, in this picture, as nucleon stars. The maximal possible mass of this type of star, 
where Eq. fl9^|) has gone to completion, has been calculated to be around 1.5 M Q [189 . Consequently, 
the collapsing core of a supernova (e.g. 1987A), if heavier than this value, should go into a black hole 
rather than forming a neutron star |172H185l l"186 . Another striking implication, pointed out by Brown 
and Bethe, would be the existence of a large number of low-mass black holes in our galaxy [172 . 
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Figure 21: Central density versus rotational frequency for several sample neutron stars. 
The stars' baryon number, A, is constant in each case. Theory predicts that the interior 
stellar density could become so great that the threshold densities of various novel phases of 
superdense matter are reached, eo = 140 MeV/fm 3 denotes the density of nuclear matter, 
f^K is the Kepler frequency, and M(0) is the stars' mass at zero rotation. (From Ref. [2].) 



5.1.3 Strange Quarks 

Already several decades ago it has been suggested that, because of the extreme densities reached in the 
cores of neutron stars, neutrons protons plus the heavier constituents may melt, creating quark matter 
being sought at the most powerful terrestrial heavy- ion colliders [171 EU 1221 122] • At present 
one does not know from experiment at what density the expected phase transition to quark matter 
occurs, and one has no conclusive guide yet from lattice QCD simulations. From simple geometrical 
considerations it follows that nuclei begin to touch each other at densities of ~ (Anrfj/3)^ 1 ~ 0.24 fm~ 3 , 
which, for a characteristic nucleon radius of r^ ~ 1 fm, is less than twice the baryon number density 
p of ordinary nuclear matter [Tj. Above this density, therefore, is appears plausible that the nuclear 
boundaries of hadrons dissolve and the originally confined quarks begin to populate free states outside 
of the hadrons. Depending on rotational frequency and stellar mass, densities as large as two to three 
times po are easily surpassed in the cores of neutron stars, as can be seen from Figs. [""""] and [221 so 
that the neutrons and protons in the centers of neutron stars may have been broken up into their 
constituent quarks by gravity [2*1 I190j . More than that, since the mass of the strange quark is so small, 
high-energetic up and down quarks are expected to readily transform to strange quarks at about the 
same density at which up and down quark deconfinement sets in "2*"H 172*] . Three flavor quark matter 
could thus exist as a permanent component of matter in the centers of neutron stars [21 [II T90J. 

In passing we mention that the existence of a mixed phase depends decisively on the unknow surface 
tension. Alford et al. studied the CFL-nuclear mixed phase |191j and found that if the surface tension is 
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Figure 22: Variation of central star density with rotational star frequency. 



above about 40 MeV/fm 2 then the mixed phase does not occur at all. This surface tension is completely 
unknown, but one mig ht expect the 'natural scale' for it to be (200 MeV) 3 = 200 MeV/fm. For what 
follows we will assume that the surface tension is such that a mixed phase of quark matter and nuclear 
matter exists above a certain density. 

As pointed out by Glendenning PUESj, in all earlier work on the quark-hadron phase transition in 
neutron star matter, assumed to be a first order transition, the possibility of reaching the lowest energy 
state by rearranging electric charge between the regions of confined hadronic matter and deconfined 
quark matter in phase equilibrium was ignored. This incorrectly yielded a description of the possible 
quark-hadron transition in neutron stars as a constant pressure one, which had the consequence of 
excluding the coexistence phase of hadrons and quarks from neutron stars. The microphysical agent 
behind this preference for charge rearrangement is the charge-symmetric nuclear force which acts to 
relieve the high isospin asymmetry of neutron star matter as soon as it is in equilibrium with quark 
matter. This introduces a net positive charge on the hadronic regions and a compensatory net negative 
charge on quark matter. (The nuclear matter and quark matter phases are thus not separately charge 
neutral as assumed before Glendenning's work.) Because of this preference for charge rearrangement 
exploited by a neutron star, the pressure in the mixed quark-hadron phase varies as the proportions of 
the phases. Varying pressure in the mixed phase is of key importance for the existence of the mixed 
phase inside neutron stars, because hydrostatic equilibrium dictates that pressure drops monotonically 
from the center toward the surface. For that reason the mixed phase is not strictly-and incorrectly- 
excluded from neutron stars . 

If the dense interior of a neutron star is converted to quark matter EU IT] I192j . it must be three- 
flavor quark matter (see Figs. ESI through l23j) since it has lower energy than two-flavor quark matter 
(see section EJ). And just as for the hyperon content of neutron stars, strangeness is not conserved on 
macroscopic time scales which allows neutron stars to convert confined hadronic matter to three-flavor 
quark matter until equilibrium brings this process to a halt. Many of the earlier investigations have 
treated neutron stars as containing only neutrons, and the quark phase as consisting of the equivalent 
number of u and d quarks. Pure neutron matter, however, is not the ground state of a neutron star, 
nor is a mixture of u and d quarks the ground state of quark matter in compact stars. In fact the latter 
constitutes a highly excited state of quark matter, which will quickly weak decay to an approximate 
equal mixture of u, d and s quarks. Several other investigations have approximated the mixed phase 
as two components which are separately charge neutral, which hides the possible quark-hadron phase 
transition in neutron star matter EH [25 1 because the deconfmement transition is shifted to densities 
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Figure 23: Composition of chemically equilibrated, stellar quark-hadron (hybrid star) matter 
as a function of baryon density. Hadronic matter is described by the relativistic Hartree 
model HV 0ES|, the bag constant is B = 250 MeV/fm 3 . (From Ref. 0.) 



hardly reached in the cores of neutron stars of average mass, M ~ 1.4 M . 

The Gibbs condition for phase equilibrium between quarks and hadrons is that the two associated, 
independent chemical potentials, \x n and \i e respectively, and the pressure in the two phases be equal, 

P H (//\ //, {0}, T) = P Q (//\ //, T) , (97) 

where /i n and fi e denote the chemical potentials of neutrons and electrons, respectively, and the sub- 
scripts H and Q refer to confined hadronic matter and deconfined quark matter. The quantity {<p} 
stands collectively for the particle fields and Fermi momenta which characterize a solution to the equa- 
tions of confined hadronic matter discussed in section [21 As known form Eq. ([10)1 . the quark chemical 
potentials are related to the baryon and charge chemical potentials in Eq. (j9Tj) as 

p u = p c =\p n -\p\ p d = f = ±n n + ±vf. (98) 

In accordance to what has bee said just above, Eq. (J9?]) is to be supplemented with the conditions 
of baryon charge conservation and electric charge conservation [2U ES]- Mathematically, the global 
conservation of baryon charge within an unknown volume, V, containing A baryons is expressed as 

P = ^ = {l-X) Ph(/A /A T) + X Pq(/A /A T) , (99) 

where x = ^q/^ denotes the volume proportion of quark matter, Vq, in the unknown volume V. By 
definition, the parameter x varies between and 1, determining how much confined hadronic matter 
exists as quark matter. The global neutrality of electric charge within the volume V is mathematically 
expressed as (2HI2S] 

= | = (1 - X) <?h(//\ /A T) + x g Q (/A /A T) + q L , (100) 

where gn = Hb ( 1'bP B an d <?q = J2flfp^ denote the net electric charge carried by hadronic and quark 
matter, respectively, and = J2l QlP L stands for the electric charge density of the leptons (see sections 
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Figure 24: Composition of chemically equilibrated, stellar quark-hadron matter as a function 
of baryon density Hyperons are artificially suppressed. Hadronic matter is described by HFV 
HIE], the bag constant is B = 150 MeV/fm 3 . (Fig. from Ref. 0.) 



El and El). One sees that for a given temperature T, Eqs. (j"""7j) through (jlOOj) serve to determine the two 
independent chemical potentials p n and p e , and the volume V for a specified volume fraction \ °f t ne 
quark phase in equilibrium with the hadronic phase. After completion, Vq is obtained as Vq = yV . 
Through Eqs. (f9Tj) to (jl00|) the chemical potentials p n and /i e obviously depend on x an d thus on 
density p, which renders all properties that depend on p n and //-from the energy density to the baryon 
and charge densities of each phase to the common pressure-density dependent, too. 

Figures 12*31 through 1231 show sample quark-lepton populations computed for representative bag con- 
stants as well as different many-body approximations employed to model confined hadronic matter 
[2*1 H9Uj . Three features emerge immediately from these populations. Firstly, one sees that the transi- 
tion from pure hadronic matter to the mixed phase occurs at rather low density of about 3 po or even 
somewhat less [2U ESI H93j . Depending on the bag constant and the underlying nuclear many-body 
approximation, threshold values even as small as about 2 po are possible. Secondly, we emphasize the 
saturation of the number of electrons as soon as quark matter appears, for then electric charge neutral- 
ity can be achieved more economically among the baryon-charge carrying particles themselves. This 
saturation is of very great importance for the possible formation of a K~ condensate in neutron stars 
|185U186tll87|ll88] . whose threshold condition is given by p e = m* K . Figure [2*01 shows that this condition 
may be fulfilled in neutron star matter, depending on the underlying many-body approximation and the 
structure of many-body background. Thirdly, the presence of quark matter enables the hadronic regions 
of the mixed phase to arrange to be more isospin symmetric, i.e. closer equality in proton and neutron 
number is achieved, than in the pure phase by transferring charge to the quark phase in equilibrium 
with it. Symmetry energy will be lowered thereby at only a small cost in rearranging the quark Fermi 
surfaces. Electrons play only a minor role when neutrality can be realized among the baryon-charge 
carrying particles. The stellar implication of this charge rearrangement is that the mixed phase region 
of the star will have positively charged regions of nuclear matter and negatively charged regions of quark 
matter. 

Because of the competition between the Coulomb and the surface energies associated with the 
positively charged regions of nuclear matter and negatively charged regions of quark matter, the mixed 
phase will develop geometrical structures, similarly as it is expected of the subnuclear liquid-gas phase 
transition |194|I195(I196| . This competition establishes the shapes, sizes and spacings of the rarer phase 
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Figure 25: Same as Fig. El but for B = 250 MeV/fm 3 . (Fig. from Ref. 0.) 



in the background of the other in order to minimize the lattice energy [2U ESI H97j . As known from the 
quark-hadron compositions shown above, the formation of quark (q) drops may set in around 3po- At 
a somewhat greater density the drops are more closely spaced and slightly larger in size. Still deeper 
in the star, the drops are no longer the energetically favored configuration but merge together to form 
quark rods of varying diameter and spacing. At still greater depth, the rods grow together to quark 
slabs. Beyond this density the forms are repeated in reverse order until at the inner edge of the mixed 
phase hadronic (h) drops of finite size but separated from each other are immersed in quark matter. 
At densities between six to ten times po, the hadronic drops have completely dissolved into pure quark 
matter [2U ESI 0]- In all cases the geometric forms lie between about 10 and 25 fm 1197] . The 
change in energy accompanied by developing such geometrical structures is likely to be very small in 
comparison with the volume energy |2"Hl2l)l HUH} 1199] and, thus, cannot much affect the global properties 
of a neutron star. Nevertheless, the geometrical structure of the mixed phase may be very important 
for transport phenomena as well as irregularities (glitches) in the timing structure of pulsar spin-down 

EHES1II!. 

We conclude this section with presenting a representative model for the equation of state of a quark- 
hybrid star, which is shown in Fig. EB1 [193 . The hadronic phase is modeled, as in Fig. ESI in the 
framework of the relativistic Hartree approximation (model HV of Refs. ffi I48| ITD~Uj ). the quark phase 
by the bag model with a bag constant of B = 150 MeV/fm 3 and a strange quark mass of 150 MeV. 
The only difference with respect to Fig. E21 is the smaller bag constant which lowers the onset of quark 
deconfinement from 3p to 2 p |193j . Up to neutron chemical potentials of p n ~ 10 3 MeV, the matter 
stays in the pure hadronic phase. The onset of quark deconfinement, which saturates the number of 
electrons (cf. Fig. I2*3*|). occurs at point a in the diagram where p e attains its maximum, /i e ~ 180 MeV. 
As remarked just above, this values corresponds to a baryon number density of about 2 p . Beyond 
this density, p e decreases toward rather small values because less and less electron are present in dense 
quark-hybrid star matter. The mixed phase region (a-b) exists, in the direction of increasing density, 
for electron chemical potential in the range 180 MeV ^ p e ^ 25 MeV, which corresponds to neutron 
chemical potential of 10 3 MeV < p n < 1.2 x 10 3 MeV. For this range the volume proportion of quark 
matter varies from < x < 1. The energy density in the mixed phase is the same linear combination 
of the two phases as the charge and baryon number |2"4*1 ESI 0] , namely 



e = (1 - x) e H (//\ /A {0}, T) + X e Q (/A p e , T) 
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Figure 26: Equation of state of neutron star matter accounting for quark deconfinement [2]. 



Most importantly, the pressure in the mixed phase region varies with density rather than being constant, 
which would be the case if the conservation of electron charge were ignored. The pure quark matter 
phase (x — 1) sets in at b where the density has grown to about 6 p . It is characterized by a relatively 
steep increase of pressure with density. Whether or not this phase exists in neutron stars constructed 
for this equation of state depends of the star's central density and thus on its mass. 



5.1.4 H-dibaryons 

A novel particle that could make its appearance in the center of a neutron star is the so-called H- 
dibaryon, a doubly strange six-quark composite with spin and isospin zero, and baryon number two 
[86J. Since its first prediction in 1977, the H-dibaryon has been the subject of many theoretical and 
experimental studies as a possible candidate for a strongly bound exotic state. In neutron stars, which 
may contain a significant fraction of A hyperons, the A's could combine to form H-dibaryons, which could 
give way to the formation of H-dibaryon matter at densities somewhere above ~ 3 eo |2UU[ I201|. 202 
depending on the in-medium properties of the H-dibaryon. For an attractive optical potential, Uu, 
of the H-dibaryon at normal nuclear density the equation of state is softened considerably, as shown 
in Fig. H-dibaryon matter could thus exist in the cores of moderately dense neutron stars. H- 
dibaryons with a vacuum mass of about 2.2 GeV and a moderately attractive potential in the medium 
of about [/h = —30 MeV, for instance, could go into a boson condensate in the cores of neutron stars 
if the limiting star mass is about that of the Hulse-Taylor pulsar PSR 1913+16, M = 1.444 M Q [2712], 
Conversely, if the medium potential were moderately repulsive, around Un = +30 MeV, the formation 
of H-dibaryons may only take place in heavier neutron stars of mass M ^ 1.6 M Q . If formed, however, 
H-matter may not remain dormant in neutron stars but, because of its instability against compression 
could trigger the conversion of neutron stars into hypothetical strange stars |201|. I203| [204 . 
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Figure 27: Equations of state of neutron star matter accounting for a H-dibaryon condensate 
202 . Uk is the optical potential of the H-dibaryon at normal nuclear density. Reprinted 
figure with permission from N. K. Glendenning and J. Schaffner-Bielich, Phys. Rev. C 58 
(1998) 1298. Copyright 1998 by the American Physical Society. 



5.2 Strange Stars 
5.2.1 General Properties 

If the strange matter hypothesis is true, a new class of compact stars called strange stars should exist. 
Possible strange star candidates as compiled in table 03 They would form a distinct and disconnected 
branch of compact stars and are not a part of the continuum of equilibrium configurations that include 
white dwarfs and neutron stars (see Fig.l28|). In principle both strange and neutron stars could coexist. 
However if strange matter is the true ground state, the galaxy is likely to be contaminated by strange 
quark nuggets which, depending on their velocities j^Gj, could convert neutron stars to strange stars 
[TU1 11521 1219j . This would mean that the objects known to astronomers as pulsars would be rotating 
strange stars (see Fig. |2*9*|) rather than rotating neutron stars. Another peculiar consequence of the 



Table 5: Possible strange star candidates. 



Compact object 


Peculiar feature 


References 


RX J1856.5-3754 


Small radius 


12051 I27T61 127171 12TTE1 127J3] 


4U 1728-34 


Small radius 


pro] 


SAX J1808.4-3658 


Small radius 


pTTj 


Her X-l 


Small radius 


|2T5] 


IE 1207.4-5209 


Peculiar timing 


|2T3] 


PSR 0943+10 


Microstorms 


eh] 


3C58 (J0205+6449) 


Low temperature 




GRO J1744-28 


X-ray burst features 


j2TB] 


SGR 0526-66, SGR 1900+14, SGR 1806-20 


X-ray burst features 


[1111111131217, 218 
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Figure 28: Mass versus radius of strange-star configurations with nuclear crusts (dashed and 
dotted curves) and gravitationally bound neutron stars and white dwarfs (solid curve). The 
strange stars carry nuclear crusts with chosen inner densities of e crus t = 4 x 10 11 g/cm 3 and 
e C rust = 10 s g/cm 3 , respectively. Crosses denote the termination points of strange-dwarfs 
sequences, whose quark matter cores have shrunk to zero. Dots refer to maximum-mass 
stars, minimum- mass stars are located at the vertical bars labeled 'b'. (Fig. from Ref. |2j.) 



hypothesis could be the existence of an entirely new class of dense white-dwarf-like strange stars, called 
strange dwarfs, plus an expansive range of planetary- like strange matter stars referred to as strange MA- 
CHOS. These objects could carry nuclear crusts that are between several hundred and several thousand 
kilometers thick ^3] • The situation is graphically illustrated in Fig. |2H1 The important astrophysical 
implication of the existence of strange MACHOS would be that they occur as natural stellar candidates 
which effectively hide baryonic matter, linking strange quark matter to the fundamental dark matter 
problem, which is currently one of the problems of greatest importance to astrophysics and cosmology. 
Observationally, the strange MACHOS could be seen by the gravitational microlensing experiments 
220., provided such objects exist abundantly enough. 

As described in section 01 strange quark matter is expected to be a color superconductor which, 
at extremely high densities, should be in the CFL phase. This phase is rigorously electrically neutral 
with no electrons required j3T]. For sufficiently large strange quark masses, however, the low density 
regime of strange quark matter is rather expected to form other condensation patterns (e.g. 2SC, CFL- 
K°, CFL-K + , CFL-7r°'~) in which electrons are present j2Hll2Zl- The presence of electrons causes the 
formation of an electric dipole layer on the surface of strange matter, as illustrated schematically in Fig. 
I3TH which enables strange quark matter stars to carry nuclear crusts jHl H3 H2] • An analytical expression 
for the electron number density and the electric field on the surface of strange quark matter can be 
derived using the Thomas-Fermi model. One finds |221j that 

9-5 xlO 35 , „ 7.2 x 10 18 Ar , , 1A „, 

n e ~ cm" 3 , E ~ - V cm" 1 , 102 

(1.2z u + 4) 3 (1.2*n + 4) 2 1 ; 

where z is a measured height above the quark surface, Z\\ = z/(10~ u cm). From Eq. (jl02|) one sees 
that very strong electric fields, on the order of ~ 10 17 V/m, may be expected near the quark surface. 
This makes it possible for a non-rotating star to support a nuclear crust with a mass up to ~ 10 5 M 
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Figure 29: Schematic illustration of the cross section of a rotating strange star carrying a 
nuclear crust. The quantities e c and e crust denote the star's central density and the density 
at the base of the crust, respectively. (Fig. from Ref. |2j.) 



[S]. The maximal possible density at the base of the crust, called the inner crust density, is determined 
by neutron drip, which occurs at about 4 x 10 11 g/cm 3 . This somewhat complicated situation of the 
structure of strange matter enveloped in a nuclear crust can be represented by a proper choice for the 
equation of state shown in Fig. [""""] The equation of state is characterized by a discontinuity in density 
between strange quark matter and nuclear crust matter across the electric dipole gap where the pressure 
of the nuclear crust at its base equals the pressure of strange matter at its surface [2~| lllj. 

Since the nuclear crust surrounding the quark matter core of a strange star is bound to the core by the 
gravitational force rather than confinement, the mass-radius relationship of strange matter stars with 
nuclear crust is qualitatively similar to the one of neutron stars and white dwarfs, which are exclusively 
bound by gravity. The strange-star sequences in Fig. [2E1 are computed for the maximal possible inner 
crust density set by neutron drip (e crus t = 4 x 10 11 g/cm 3 ) as well as for a considerably smaller sample 
density of 10 8 g/cm 3 . Of course there are other possible sequences of strange stars with any smaller value 
of inner crust density [21 ""Q1 [E] • From the maximum-mass star (solid dots), the central star density 
decreases monotonically through the sequence in each case. The fact that strange stars with nuclear 
crusts possess smaller radii than neutron stars leads to smaller Kepler (mass shedding) periods, Pk, for 
strange stars. This is indicated by the classical expression Pk oc \JbFJm and has its correspondence in 
the general relativistic expression for P K derived in Eq. (J"""">J). Since the qualitative dependence of Pk 
on mass and radius remains valid [222], one finds that the complete sequence of strange stars, and not 
just those close to the mass peak as is the case for neutron stars, can sustain extremely rapid rotation. 
In particular, a strange star with a typical pulsar mass of ~ 1.4 M & can rotate at Kepler periods as 
small as Pk ~ 0.5 ms, depending on crust thickness and the modeling of strange quark matter This 
is to be compared with the larger limiting value of Pk ~ 1 ms obtained for neutron stars of the same 
mass [2 j . Exceptions to this, however, are possible if the nuclear matter exhibits a very strong softening 
at intermediate densities |177j . 

Of considerable relevance for the viability of the strange matter hypothesis is the question of whether 
strange stars can exhibit glitches in rotation frequency. Pulsar glitches are sudden changes in the 
rotational frequency of a rotating neutron star which otherwise decreases very slowly with time due to 
the loss of rotational energy through the emission of electromagnetic dipole radiation and an electron- 
positron wind. They occur in various pulsars at intervals of days to months or years, and in some 
pulsars are small (Crab), and in others large (Vela) and infrequent (AQ/Q ~ 1CT 8 — 1CT 6 , respectively). 
Glitches have been attributed to several causes related to the assumed structure of neutron stars. One 
such is the crust quake in which an oblate solid crust in its present shape slowly comes out of equilibrium 
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Figure 30: Surface region of strange quark matter. The region between R < r < R crns t 
is filled with electrons that are bound to strange matter but extend beyond its surface, R, 
leading to a deficit of electrons in the range R m < r < R and therefore a net positive charge 
in this region. The associated electric field, E ~ 10 17 V/cm, is sufficiently strong to avoid 
contact between atomic matter and strange matter, enabling strange matter to be enveloped 
by ordinary atomic matter. (Fig. from Ref. |2j.) 



with the forces acting on it as the period of rotation changes, and fractures when the built up stress 
exceeds the sheer strength of the crust material 223 , 224]. The stellar frequency and rate of change of 
frequency, Q and Cl respectively, slowly heal to the trend preceding the glitch as the coupling between 
crust and core reestablishes their corotation. The compatibility of pulsar glitches with the strange 
matter hypothesis will have a decisive impact on the question of whether or not strange matter is the 
true ground state of strongly interacting matter [TT| !22o| . From the study performed in Refs. j^, 226 
it is known that the ratio of the crustal moment of inertia to the total moment of inertia, / C rust/-f total, 
varies between 10 -3 and ~ 10~ 5 . If the angular momentum of the pulsar is conserved in a stellar quake 
then the relative frequency change and moment of inertia change are equal, and one arrives for the 
change of the star's frequency at [11 J 

^ = ^>^ = /^~(10- 5 -10- 3 )/, with 0</<l. (103) 

Here Iq denotes the moment of inertia of that part of the star whose frequency is changed in the quake. 
It might be that of the crust only, or some fraction, or all of the star. The factor / in Eq. ()103|) represents 
the fraction of the crustal moment of inertia that is altered in the quake, i.e., / = |AJ|/J crust . Since 
the observed glitches have relative frequency changes AQ/Q = (10~ 9 — 10~ 6 ), a change in the crustal 
moment of inertia of / & 0.1 would cause a giant glitch even in the least favorable case [IT] . Moreover 
it turns out that the observed range of the fractional change in the spin-down rate, Q, is consistent with 
the crust having the small moment of inertia calculated and the quake involving only a small fraction 
/ of that, just as in Eq. ()103|) . For this purpose we write [TT] 

where use of Eq. (J103J) has been made. Equation ()104|) yields a small / value, i.e., / < (10~ 4 to 10 _1 ), in 
agreement with / £ 10 _1 established just above. Here measured values of the ratio (AO/0)/(Af2/0) ~ 
10~ 6 to 10~ 4 for the Crab and Vela pulsars, respectively, have been used. 
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Figure 31: Equation of state of strange quark matter surrounded by a nuclear crust [TT] . 
The maximal possible nuclear matter density is determined by neutron drip which occurs 
at e C mst — 0.24 MeV/fm 3 (4.3 x 10 11 g/cm 3 ). Any nuclear density that is smaller than 
neutron drip is possible. As an example, we show the equation of state for a sample density 
of e crust = 10" 4 MeV/fm 3 (10 8 g/cm 3 ). 



5.2.2 SAX J1808. 4-3658 

As discussed just above, strange stars are self-bound objects at zero external pressure, which would exist 
stably even if gravity were switched off. The latter makes such objects even more compact. The radii 
of compact quark stars are thus expected to be smaller than those of neutron stars of comparable mass, 
which are bound solely by gravity. Figures CE1 arid E21 show that this difference in radius may be 
three to four kilometers for stars of canonical mass, ~ 1.4 M , and even bigger than that for less massive 
objects. One neutron star that may have such an unusually small radius is the transient X-ray burst 
source SAX J1808. 4-3658 |211j . which was discovered in September 1996 by the BeppoSAX satellite. 
Two bright type-I X-ray bursts were detected, each one lasting less than 30 seconds. An analysis of 
the bursts in SAX J 1808. 4-3658 indicates that it is 4 kpc distant and has a peak X-ray luminosity of 
6 x 10 36 erg/s in its bright state, and a X-ray luminosity lower than 10 35 erg/s in quiescence [227J. 
Coherent pulsations at a period of 2.49 milliseconds were also discovered [228J. The binary nature of 
SAX J1808. 4-3658 was firmly established with the detection of a 2 hour orbital period |229j as well as 
with the optical identification of the companion star. SAX J1808. 4-3658 is the first pulsar that shows 
both coherent pulsations in its persistent emission and X-ray bursts. 

Li et al. j211j extracted a mass-radius relationship for the compact star in SAX J1808. 4-3658 from the 
following two requirements on the geometry of the stellar binary system (see also section I7~2l where mass 
accretion onto a compact star is discussed in more detail). Firstly, detection of X-ray pulsations requires 
that the inner radius, r m , of the accretion disk should be larger than the stellar radius R. Secondly, the 
inner radius must be smaller than the disk's co-rotation radius, r c . Otherwise accretion will be inhibited 
by a centrifugal barrier. From these two conditions one finds that r m < r c = (MP 2 /(4ff 2 )) 1/3 , with M 
and P the star's mass and rotational period, respectively. Expressing the location of the inner disk in 
terms of the Alfven radius, ta, which is that distance from the neutron star at which the accreting matter 
is pulled by the magnetic fields to the poles of the star which happens when the kinetic energy density 
is comparable to the magnetic energy density |230j . one arrives at r m = ^r^ = £(B 2 R 6 / M{2M) 1 / 2 ) 2 ^ . 
The symbols B and M are the surface magnetic field and the mass accretion rate of the rotating neutron 
star, respectively. The quantity £ ~ 1 is a parameter which depends very weakly on the accretion rate 
[231 . Denoting the minimum and maximum accretion rates of SAX J1808. 4-3658 as M m ; n and M max , 
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Figure 32: Comparison of the mass-radius relationship of SAX J1808. 4-3658 determined from 
RXTE observations (i.e. Eq. ([106)) ) with theoretical models of neutron stars and strange stars 
|211j . The solid curves represents theoretical mass-radius relationships for neutron stars and 
strange stars. Reprinted figure with permission from X. D. Li et ai, Phys. Rev. Lett. 83 
(1999) 3776. Copyright 1999 by the American Physical Society. 



the conditions discussed above can be expressed as 

r < r m (M max ) < r m (M min ) < r c . (105) 

To connect this relation to the observed data, let us assume that the mass accretion rate is proportional 
to the X-ray flux F observed with the Rossi X-ray Timing Explorer (RXTE). This is supported by 
the fact that the X-ray spectrum of SAX J1808. 4-3658 was remarkably stable and that there was 
only a slight increase in the pulse amplitude when the X-ray luminosity varied by a factor of ~ 100 
during the 1998 April/May outburst [232, 1233"! I234j . From this, an upper limit for the stellar radius of 
R < (F min / F maiX ) 2 / 7 r c was derived in |211j . This relation can be conveniently written as 

where -F max and F min denote the X-ray fluxes measured during high and low X-ray emission states, 
respectively. Adopting a flux ratio of F max / F min ~ 100, Eq. ()106|) constrains the mass-radius values 
of SAX J1808. 4-3658 to values that lie between the dashed curves in Fig. E21 which suggests that 
SAX J1808. 4-3658 could be a strange star. The dashed line labeled R = R s (= 2MR) denotes the 
Schwarzschild limit on the radius of a compact object set by gravitational collapse to a black hole. 
The curves labeled BBBl, BBB2, Hyp, and K~ denoted the mass-radius relationships of conventional 
neutron stars computed for different models for the equation of state. Constraints on the mass-radius 
relationship which, if robust, are better described in terms of strange stars than neutron stars have also 
been obtained for the compact star in the X-ray source 4U 1728-34 (M < 1.0 M & , R < 9 km) [235], for 
the isolated neutron star RX J1856. 5-3754 |2()5j . which will be discussed in the next section, as well as 
for the X-ray pulsar Her X-l (M = 1.1 - 1.8 M , R = 6.0 - 7.7 km) (213 • (See also table 0) 

In passing, we mention the recent discovery of significant absorption lines in the spectra of 28 bursts 
of the low-mass X-ray binary EXO 0748-1676 23^1 . These lines have been identified with iron and 
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Figure 33: The critical temperature for hydrogen (H) and iron (Fe) as a function of magnetic 
field, B 205J. A phase transition to a solid condensate is possible in the shaded region for Fe 
and in the cross-hatched region for H. The positions of five cool, isolated neutron stars (see 
table |7j) is shown. The horizontal line is drawn in correspondence to the color temperature 
of RX J1856. 5-3754. Reprinted figure with permission from R. Turolla, S. Zane, and J. J. 
Drake, Astrophys. J. 603 (2004) 265. Copyright 2004 by the Astrophysical Journal. 

oxygen transitions, all with a gravitational redshift of z = 0.35. As shown in Ref. |236j . for a stellar 
mass range of M ~ 1.3 — 2.0 M such a z value is completely consistent with conventional neutron star 
models, made of normal nuclear matter, and excludes even some models in which neutron stars are 
made of more exotic matter. 

5.2.3 RX J 1856. 5-3754 

Aside from X-ray emission from neutron stars in binaries, discussed in the previous section, the thermal 
emission from the surface of an isolated neutron star (INS) is of key importance for the determination 
of the mass and radius of a neutron star, too. ROSAT was the first satellite with a sufficient sensitivity 
in the X-ray band to perform a systematic search and study of such objects. The nearby neutron 
star RX J1856. 5-3754, discovered in 1992 |237j . is the brightest INS in X-rays. It does not show any 
signs of activity such as variability or pulsation. Since its discovery RX J1856. 5-3754 was studied in 
great detail in X-rays, UV and the optical band using and variety of different astrophysical observa- 
tories (ROSAT, EUVE, ASCA, HST, Chandra, XMM Newton, VLA). Detailed Chandra observations 
of RX J1856. 5-3754 have shown that this neutron star has a featureless thermal spectrum for which 
a simple blackbody distribution seems to provide a better fit to X-ray data than more sophisticated 
atmospheric models |238| [239 . Several sets of observations taken with ROSAT, EUVE, NTT, Keck, 
HST and Chandra have yielded a proper motion of 330 mas/yr and parallax 8.5 mas, which corre- 
sponds to a transverse stellar velocity of about 200 km/s and a distance from Earth of D = 120 pc 
[240, 1241] . These measurements could indicate that the star originated ~ 10 6 years ago in the Upper 
Scorpio association at about the same time that a supernova ejected the runaway O star ( Ophiuchus, 
suggesting that RX J1856. 5-3754 and £ Ophiuchus could have been members of a binary system and 
that RX J1856. 5-3754 may be the remnant of the star that exploded. The apparent (redshifted) radius 
of RX J1856.5-3754 can be calculated from [H| 



IT = 4.25 fX 1/2 fE 1/2 V Dioo TC/60 eV)- 2 km. 



(107) 
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The true (coordinate) radius of the star, R, is connected to the apparent radius through 

R = R°°(l -2M/R) l/2 . (108) 

The quantity J'a in Eq. ()107j) denotes the fraction of the stellar surface responsible for the observed 
emission, fg is the ratio of emitted power to blackbody power, 7 stands for the spectral hardening, 
D wo = D I (100 pc) is the star's distance, and T bb = T{£(1 — 2M/i?) -1 / 2 its blackbody temperature. If 
the emission comes from the entire stellar surface (Ja = 1) and the temperature distribution is uniform 
(7 = 1), Eq. ()107|) simplifies to 

R°° = A.25f E 1/2 D 100 (T bb /60 eV)~ 2 km . (109) 

Based on this relation Drake et al. 206] reported an apparent radius for RX J1856. 5-3754 of R°° / D wo = 
4.12 ± 0.68 km, where the quoted uncertainty represents the combined temperature determination 
uncertainty (±1 eV) and the dominant absolute effective area uncertainty of Chandra's Low Energy 
Transmission Grating (LETG) and High Resolution Camera spectroscopic plate detector array (HRC-S) 
combination (±15%). Depending on the star's distance, which still appears to be an open issue [242J, 
this implies a rather small apparent stellar radius in the range between about R°° = 4 km to 8 km. The 
latest distance measurement of 175 pc [242 combined with specific assumptions about the star's surface 
composition (see below) could favor _R°° ~ 8 km |243j . which corresponds to a true radius of about 
6 km. As can be seem from Figs. EU and EH such a value would be too small for conventional neutron 
star equations of state 1206] which predict 12 km £ R°° £ 17 km (3j. Proposed explanations for such 
a small radius include a reduced X-ray emitting region (such as a heated polar cap), or the presence 
of a more compact object such as a strange quark star |207t I208| I209j . An alternative possibility, 
suggested in |205j . could be cold neutron stars (T £ 10 6 K) endowed with rather strong magnetic fields 
(B 10 13 G) and with metal dominated outer layers. Depending on T and B such stars may undergo a 
phase transition to a solid condensate in the outermost layers (Fig. I33j) . resulting in an X-ray spectrum 
that is featureless as observed for RX J1856. 5-3754. The observed UV-optical enhanced emission could 
be explained by the presence of a gaseous, thin hydrogen shell on top of the iron condensate, where 
the optical flux is reprocessed. This model predicts a value for the apparent radius of RX J1856. 5-3754 
in the range of R°° ~ 8 — 12 km, depending on whether one assumes uniform or meridional variation 
temperature distributions on the star's surface. 

For a canonical neutron star mass of 1.4 M Q , apparent radii in the range R°° ~ 10 — 12 km (true 
radii of R ~ 7 — 9 km for a 1.4 M neutron star) can be accommodated by conventional models for the 
equations of state of superdense matter which exhibit a soft behavior at high densities |205j . This is 
different for R°° £ 8 km which may indicate that such stars are made of self-bound hadronic matter, of 
which strange quark matter may be the most plausible state of matter. Sequences of such stars are shown 
in Figs. El and E21 While a strange star may be an intriguing conceivable option, present observations of 
RX J1856. 5-3754 do not necessarily demand this interpretation and more conventional interpretations 
involving conventional neutron stars are certainly possible |205t I244j . A conventional interpretation is 
also supported by models which fit the X-ray and optical data of RX J1856. 5-3754 with a two-component 
blackbody model [233] . The latter are best fitted with T{g x ~ 63.5 eV and ~ 4.4(D/120 pc) km 

for the hot X-ray emitting region, and t < 33 eV and R^ opt > 17(Z?/120 pc) km for the rest of 
the neutron star surface responsible for the optical flux |239j . 

5.2.4 The Neutron Star in 3C58 

Murray et al. discovered the 65 ms pulsars PSR 0205+6449 in the supernova remnant 3C58 located 
in the constellation Cassiopeia |245j . Historical evidence strongly suggests an association of the rem- 
nant with supernova SN 1181, which went off in 1181 AD. This renders PSR 0205+6449 younger than 
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B0531+21, the pulsar in the Crab nebula, born in a supernova explosion recorded by Chinese as- 
tronomers in 1045 AD. Using data provided by the Chandra X-ray observatory, Slane et al. inferred an 
upper limit on the effective surface temperature of PSR 0205+6449of only T°° < 1.08 x 10 6 K which falls 
below predictions from standard cooling models |215j . Based on its low temperature, it was suggested 
that PSR 0205+6449 may be a strange quark star rather than a neutron star |215j . We will discuss 
this issue in greater detail in section |3 

5.2.5 X-ray, Gamma-Ray Burst and SGR Associations 

X-ray bursts are sudden increases in the X-ray flux of X-ray sources, with rise times of < Is, and 
subsequent decay with characteristic times ranging from 10s to a few minutes. They are classified into 
two classes: type-I and type-II X-ray bursts. Type-I X-ray bursts are characterized by the relatively 
long burst intervals (hours to days) and the significant spectral softening during the burst decay as 
compared to their type-II counterparts. The differences between type-I(II) bursts and the HXRBs is 
discussed by [246 . 

Gamma ray bursts (GRBs) are most intense transient gamma ray events in the sky when they are 
on but the nature of gamma ray bursts is still a mystery since their discovery [247J. The isotropic and 
inhomogeneous distribution of GRBs detected with BATSE and the identification of counterparts in 
radio, optical, and X-ray by BeppoSAX and other ground based telescopes supports that they are located 
at cosmological distances, which make them the most energetic events ever known [2482 1249j . Gamma 
ray bursts at truly cosmological distances could be due to collisions of two neutron or two strange stars 
in binary systems 250 , depending on the true ground state of strongly interacting matter, and/or may 
also involve black holes |251j . Binary neutron star or strange star collisions could release ~ 10 52 ergs in 
the form of gamma rays over a time period of about 0.2 s. The central engine that powers gamma ray 
bursts should be capable of releasing a total energy of ~ 10 53 erg, which may or may not be possible 
in stellar collisions. This may be different for the conversion of a neutron star to a strange star 252J. 
Depending on the model for the equation of state the total energy given off in such a conversion is 
(1 - 4) x 10 53 erg |252| . 

Bare strange stars have also been associated with soft gamma repeaters (SGRs). Unlike gamma ray 
bursts, which emit large amounts of high energy gamma radiation, SGRs have a larger proportion of 
lower-energy X-ray radiation. Also, in contrast to gamma ray bursts, which can rumble on for many 
minutes, SGRs pop off in a fraction of a second. Conventional models associate SGRs with young 
neutron stars energizing a large cloud of gas cast off in a supernova explosion, or X-ray binary stars 
that accrete matter at irregular intervals and emit gamma rays when the accreted matter hits the 
surface. Most strikingly, the intensity of the outbursts is between 10 3 and 10 4 times the Eddington 
limit, ^Edd- The latter is defined as the critical luminosity at which photon radiation pressure from the 
surface of a star of mass M equals gravity, 

L Edd = 1.3 x 10 38 (M/M Q ) erg s" 1 . (110) 

The Eddington limit does not apply to bare strange stars since quark matter is held in place by the 
strong interaction and not gravity (see section lo.2.7|) . For that reason the bursting activity of SGRs 
can be comfortably explained by fast heating of the surfaces of bare quark stars up to temperatures of 
~ (1 — 2) x 10 9 K and its subsequent thermal emission jlUl^ 1112] . The fast heating mechanism of SGRs 
may be either impacts of comets onto bare strange stars |1011 11131 1114j or fast decay of superstrong 
magnetic fields |253l 1115] . 

A very high-luminosity flare took place in the Large Magellanic Cloud (LMC), some 55 kpc away, 
on 5 March 1979. Another giant flare was observed on 27 August 1998 from SGR 1900+14. The 
inferred peak luminosities for both events is ~ 10 7 (i.e. ~ 10 45 erg/s) times the Eddington limit for 
a solar mass object, and the rise time is very much smaller than the time needed to drop ~ 10 25 g 
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Figure 34: Temperature dependence of the critical angular velocity Qc of rotating neutron 
stars. The left figure shows the gravitational radiation driven /-mode instability suppressed 
by shear and bulk viscosity (Fig. from Ref. Right figure: comparison of /-mode insta- 
bility with r-mode instability. (Data from Refs. |2551 [256 .) 



(about 10~ 8 M Q ) of normal material onto a neutron star. Alcock et al. jH] suggested a detailed model 
for the 5 March 1979 event burst which involves the particular properties of strange matter (see also 
Horvath et al. |101[ 111 Jj ). The model assumes that a lump of strange matter of ~ 1O~ 8 M fell onto 
a rotating strange star. Since the lump is entirely made up of self-bound high-density matter, there 
would be only little tidal distortion of the lump, and so the duration of the impact can be very short, 
around ~ 10~ 6 s, which would explain the observed rapid onset of the gamma ray flash. The light curves 
expected for such giant bursts |112t 11141 11161 H17j should posses characteristic features that are well 
within the capabilities of ESA's INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL 
|254j ) launched by the European Space Agency in October of 2002. 

5.2.6 Rotational Instabilities 

An absolute limit on rapid rotation is set by the onset of mass shedding from the equator of a rotating 
star. However, rotational instabilities in rotating stars, known as gravitational radiation driven insta- 
bilities, set a more stringent limit on rapid stellar rotation than mass shedding. These instabilities 
originates from counter-rotating surface vibrational modes which at sufficiently high rotational star 
frequencies are dragged forward. In this case gravitational radiation, which inevitably accompanies the 
aspherical transport of matter, does not damp the instability modes but rather drives them. Viscosity 
plays the important role of damping these instabilities at a sufficiently reduced rotational frequency 
such that the viscous damping rate and power in gravity waves are comparable. The most critical 
instability modes that are driven unstable by gravitational radiation are /-modes and r-modes. Figure 
13*41 shows the stable neutron star frequencies if only /-modes were operative in neutron star. One sees 
that hot as well as cold neutron stars can rotate at frequencies close to mass shedding, because of the 
large contributions of shear and bulk viscosity, respectively, for this temperature regime. The more 
recently discovered r-mode instability may change the picture completely, as can be seen from Fig. 
These modes are driven unstable by gravitational radiation over a considerably wider range of angular 
velocities than the /-modes (cf. dashed curve labeled (m = 2) r-mode instability). In stars with cores 
cooler than ~ 10 9 K, on the other hand, the r-mode instability may be completely suppressed by viscous 
phenomena so that stable rotation would be limited by the /-mode instability again [256J. 

Figures ESI and ESI are the counterparts to Fig. but calculated for strange stars made of CFL and 
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Figure 35: Critical rotation frequen- 
cies versus stellar temperature for CFL 
strange stars |258j . Reprinted figure with 
permission from J. Madsen, Phys. Rev. 
Lett. 85 (2000) 10. Copyright 2000 by 
the American Physical Society. 



Figure 36: Same as Fig. ESJ but for 2SC 
quark stars |258j . Reprinted figure with 
permission from J. Madsen, Phys. Rev. 
Lett. 85 (2000) 10. Copyright 2000 by 
the American Physical Society. 



2SC quark matter, respectively |257U258j . The r-mode instability seems to rule out that pulsars are CFL 
strange stars, if the characteristic time scale for viscous damping of r-modes are exponentially increased 
by factors of ~ A/T as calculated in 257]. An energy gap as small as A = 1 MeV was assumed. For 
much larger gaps of A ~ 100 MeV, as expected for color superconducting quark matter (see section 
13. 2|) . the entire diagram would be r-mode unstable. The full curve in Fig. EH is calculated for a strange 
quark mass of m s = 200 MeV, the dotted curve for m s = 100 MeV. The box marks the positions of most 
low mass X-ray binaries (LMXBs) |259j . and the crosses denote the most rapidly rotating millisecond 
pulsars known. All strange stars above the curves would spin down on a time scale of hours due to 
the r-mode instability, in complete contradiction to the observation of millisecond pulsars and LMXBs, 
which would rule out CFL quark matter in strange stars (see, however, Ref. |260j ). Figure 13*61 shows the 
critical rotation frequencies of quark stars as a function of internal stellar temperature for 2SC quark 
stars. For such quark stars the situation is less conclusive. Rapid spin-down, driven by the r-mode 
gravitational radiation instability, would happen for stars above the curves. 



5.2.7 Surface Properties of Strange Stars 

Strange quark matter with a density of about two times the density of nuclear matter may exist up to 
the surface of a strange star jS] . Such a bare strange star differs qualitatively from a neutron star which 
has a density at the surface of about 0.1 to 1 g/cm 3 . As know from section 15.2. II the thickness of the 
quark surface would be just ~ 1 fm, the length scale of the strong interaction. Electrons are held to 
quark matter electrostatically, and the thickness of the electron surface is several hundred fermis. Since 
neither component, electrons and quark matter, is held in place gravitationally, the Eddington limit 
(Eq. (jllOj) ) to the luminosity that a static surface may emit does not apply, and thus the object may 
have photon luminosities much greater than 10 38 erg/s. It was shown by Usov |116j that this value may 
be exceeded by many orders of magnitude by the luminosity of e + e~ pairs produced by the Coulomb 
barrier at the surface of a hot strange star. For a surface temperature of ~ 10 11 K, the luminosity in 
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the outflowing pair plasma was calculated to be as high as ~ 3 x 10 51 erg/s. Such an effect may be a 
good observational signature of bare strange stars |112t 11141 Ultit I117j . If the strange star is enveloped 
by a nuclear crust however, which is gravitationally bound to the strange star, the surface made up of 
ordinary atomic matter would be subject to the Eddington limit. Hence the photon emissivity of such 
a strange star would be the same as for an ordinary neutron star. If quark matter at the stellar surface 
is in the CFL phase the process of e + e~ pair creation at the stellar quark matter surface may be turned 
off, since cold CFL quark matter is electrically neutral so that no electrons are required and none are 
admitted inside CFL quark matter [HI] . This may be different for the early stages of a hot CFL quark 
star pHTj . 

5.3 Strange Dwarfs 

The strange white dwarfs constitute the strange counterparts of ordinary white dwarfs. They consist 
of a strange quark matter core in the star's center which is enveloped by a nuclear curst made up 
of ordinary atomic matter jS]. The crust is suspended out of contact with the quark core due to the 
existence of an electric dipole layer on the core's surface [El [12]; which insulates the crust from conversion 
to quark matter. Even so, the maximum density of the crust is strictly limited by the neutron drip 
density (e drip = 4 x 10 11 g/cm 3 ), at which neutrons begin to drip out of the nuclei and would gravitate 
into the core where they would be dissolved into strange matter. 

Strange white dwarfs comprise a largely unexplored consequence of the strange matter hypothesis. 
Depending on the amount of crust mass, their properties may differ considerably from those of ordinary 
white dwarfs. For instance, it is well known that the maximum density attained in the limiting-mass 
white dwarf is about e W( j = 10 9 g/cm 3 262"l l2fci3| . Above this density, the electron pressure is insufficient 
to support the star, and there are no stable equilibrium configurations until densities of the order of 
nuclear density (e ^ 10 14 g/cm 3 ) are reached, which are neutron or strange stars. One class of strange 
dwarfs can be envisioned as consisting of a core of strange matter enveloped within what would otherwise 
be an ordinary white dwarf. They would be practically indistinguishable from ordinary white dwarfs. 
Of greater interest is the possible existence of a new class of white dwarfs that contain nuclear material 
up to the neutron drip density, which would not exist without the stabilizing influence of the strange 
quark core ^3] • The density at the inner edge of the nuclear crust carried by these strange dwarfs could 
fall in the range of e w( j < e crust < edrip- The maximum inner crust density therefore could be about 400 
times the central density of the limiting-mass white dwarf and 4 x 10 4 times that of the typical 0.6 M & 
white dwarf. An investigation of the stability of such very dense dwarf configurations to acoustical 
(radial) vibrations, which will be discussed below, reveals stability over an expansive mass range from 
~ 10~ 3 M© to slightly more than 1 M & . This is the same range as ordinary white dwarfs except that the 
lower mass limit is smaller by a factor of ~ 10~ 2 . This is because of the influence of the strange core, 
to which the entire class owes its stability. Whether or not a star is stable against radial oscillations is 
determined by a stability analysis against radial oscillations l 26^ \2bb\ . The adiabatic motion of a star 
in its nth normal eigenmode (n = is the fundamental mode) is expressed in terms of an amplitude 
f„(r) given by 



where 5r(r,t) denotes small Lagrangian perturbations in r. The quantity u n is the star's oscillation 
frequency, which we want to compute. The eigenequation for £ n (r), which governs the normal modes, 
is of Sturm-Liouville type, 



5r(r, t) 



r 



-2 



(111) 




(112) 
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Figure 37: Oscillation frequencies of the 
lowest four (n = 0,1,2,3) normal radial 
modes of strange stars with e crust = e^ip 
function of central star density . The 
cross at 'd' refers to the termination point 
of the strange dwarf sequence shown in 
Fig. El (Fig. from Ref. 
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Figure 38: Same as Fig. EH but for strange 
dwarfs in the vicinity of the termination 
point, 'd', of the sequence. The labels 'b' 
and 'c' refer to the lightest and heaviest 
star, respectively, marked in Fig. 1281 (Fig. 
from Ref. [12].) 



which implies that the eigenvalues uu 2 are all real and form an infinite, discrete sequence Uq < uo\ < 
u)\ < .... Another consequence is that the eigenfunction £ n corresponding to u 2 has n nodes in the 
radial interval < r < R. Hence, the eigenfunction £ is free of nodes in this interval. The functions 
Il(r), Q(r), and W(r) are expressed in terms of the equilibrium configurations of the star by 



n = e 2A+H r 2 rp, 

Q = -4e 2A+ ^r- 3 ^-8ne^ A+2 ^r- 2 P(e + P) + e 2A+ ^r- 2 (e + P)- 1 (^-\ 

dr \dr J 

W = e 6A+2 *r- 2 (e + P) . 



The quantities e and P denote the energy density and the pressure of the stellar matter. The pressure 
gradient, dP/dr, is obtained from the Tolman-Oppenheimer-Volkoff Eq. (J75j) . The symbol T denotes 
the adiabatic index at constant entropy, s, given by 

_ d\n P(p, s) _ (6 + P) dP(e,8) 

dlnp I' <)< 1 '' 

which varies throughout the star's interior. Solving (|112j) subject to the boundary conditions £ n oc r 3 at 
the star's origin and d£ n /dr = at the star's surface leads to the ordered frequency spectrum u 2 < uj 2 +1 
(n = 0,1,2,.. .) of the normal radial stellar modes. If any of these is negative for a particular star, the 
frequency is imaginary to which there corresponds an exponentially growing amplitude of oscillation. 
Such stars are unstable. Figures IBTI and I3~B1 show the solutions to Eq. ()112|) for the strange star sequence 
in Fig. 1251 whose inner crust density is equal to neutron drip. It follows that all stars between 'c' and 
'd' are unstable against radial pulsations since for them uo 2 < 0. This is different for the stars to the 
left of 'c' whose eigenfrequencies are all positive [T2] . 

At present there is neither a well studied model for the formation of hypothetical strange dwarfs, 
nor exists a study that determines their abundance in the universe. One possible scenario would be the 
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formation of strange dwarfs from main sequence progenitors that have been contaminated with strange 
nuggets over their lifetimes. The capture of strange matter nuggets by main sequence stars is probably 
an inevitable consequence if strange matter were more stable than hadronic matter |118j because then 
the Galaxy would be filled with a flux of strange nuggets which would be acquired by every object they 
come into contact with, like planets, white dwarfs, neutron stars and main sequence stars. Naturally, 
due to the large radii of the latter, they arise as ideal large-surface long integration time detectors for 
the strange matter flux |257j . Nuggets that are accreted onto neutron stars and white dwarfs, however, 
never reach their centers, where the gravitational potential is largest, because they are stopped in the 
lattice close to the surface due to the large structural energy density there. This prevents such stars 
from building up a core of strange matter. The situation is different for main sequence stars which 
are diffuse in comparison with neutron stars and white dwarfs. In this case the accreted nuggets may 
gravitate to the star's core, accumulate there and form a strange matter core that grows with time until 
the star's demise as a main sequence star ^21 IT3j . An upper limit on the baryon number of strange 
matter that may accumulate in a main sequence star is given by fUJ I118j 

A = 1.6 x 10 47 (M/M & )- - 15 v^ p 2i (l + O.164t; 2 2 5O (M/M )- a25 ) , (115) 

with t> 250 = foo/(250 km/s) and p 2 4 = Poo/(10~ 24 g/cm 3 ). The quantities and are the nugget 
speed and contribution to the density of the galactic halo far from the star, respectively. If one assumes 
that all the dark matter in the halo of our galaxy would consist of strange nuggets (certainly a crude 
overestimate) then v^o P24 ~ 1> and Eq. (|115|) leads for typical progenitor star masses of ~ 1 M Q to the 
upper limit of A ~ 10 48 . The mass and radius of such a strange core are ~ 2 x 10 21 kg (~ 1O _9 M 
or about 10 -3 times the mass of the Earth) and ~ 11 m. These values follows from the approximate 
relations R = 4.35 x 10 26 (M/kg) and R/km = 1.12 x KT 18 ^/ 3 for strange matter with B 1 ^ = 145 
MeV. 

Another plausible mechanism has to do with primordial strange matter bodies. Such bodies of masses 
between 10~ 2 and about 1 M & may have been formed in the early universe and survived to the present 
epoch 266J. Such objects will occasionally be captured by a main sequence star and form a significant 
core in a single and singular event. The core's baryon number, however, cannot be significantly larger 
than ~ 5 x 10 31 (M/M Q ) _L8 where M is the star's mass. Otherwise a main sequence star is not capable 
of capturing the strange matter core |I0]. Finally we mention the possibility that in the very early 
evolution of the universe primordial lumps 1119] of hot strange matter will have evaporated nucleons 
which are plausibly gravitationally bound to the lump. The evaporation will continue until the quark 
matter has cooled sufficiently Depending on the original baryon number of the quark lump, a strange 
star or dwarf, both with nuclear crusts could have been formed. 

For many years only rather vague tests of the theoretical mass-radius relationship of white dwarfs 
were possible. Recently, however, the quality and quantity of observational data on the mass-radius 
relation of white dwarfs has been reanalyzed and profoundly improved by the availability of Hipparcos 
parallax measurements of several white dwarfs |267j . In that work Hipparcos parallaxes were used to 
deduce luminosity radii for 10 white dwarfs in visual binaries of common proper-motion systems as well 
a 11 field white dwarfs. Since that time complementary follow-up HST observations have been made 
[268, 1269j to better determine the spectroscopy for Procyon B and pulsation of G226-29. Procyon B 
at first appeared as a rather compact star which, however, was later confirmed to lie on the normal 
mass-radius relation of white dwarfs. Stars like Sirius B and 40 Erin B, fall nicely on the expected 
mass-radius relation too. Several other stars of this sample (e.g. GD 140, EG 21, EG 50, G181-B5B, 
GD 279, WD2007-303, G238-44) however appear to be unusually compact and thus could be strange 
dwarf candidates |270j . 
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Figure 39: A comparison of v e count rates expected in SuperK from a proto- neutron star 
containing either np or np + Q matter |275| . Reprinted figure with permission from J. A. 
Pons et al, Phys. Rev. Lett. 86 (2001) 5223. Copyright 2001 by the American Physical 
Society. 



6 Neutrino Emission and Stellar Cooling 

The detection of thermal photons from the surfaces of neutron stars via X-ray observatories serves as 
the principal window on the properties of such objects. The surface temperatures of neutron stars are 
derivable from the measured photon flux and spectrum. The predominant cooling mechanism of hot 
(temperatures of several ~ 10 10 K) newly formed neutron stars immediately after formation is neutrino 
emission. Immediately after the birth of a (proto) neutron star in a supernova explosion, neutrinos are 
trapped inside the star because their mean free paths are shorter than the stellar radius. About 10 
seconds after birth most of the neutrinos have left the star by diffusion. As shown in |271[ 12721 1273j the 
possible presence of quark matter in the core of a neutron star could modify the diffusion rate slightly. 
Depending on distance and stellar mass the neutrino bursts will be detectable by terrestrial neutrino 
detectors such as SuperK, 1MB, Kamioka, SNO, and UNO |274j . The v e count rates expected for SuperK 
from proto-neutrons star containing either only nucleons (np) or nucleons plus quark matter (np + Q) 
are shown in Fig. EH1 where the left panel shows times less than 10 s, while the right panel shows times 
greater than 10 s. Fig. 1401 shows the total neutrino luminosity for proto- neutron stars containing quark 
matter in their centers. The shaded bands illustrate the limiting luminosities corresponding to a count 
rate of 0.2 Hz, assuming a supernova distance of 50 kpc for 1MB and Kamioka, and 8.5 kpc for SNO and 
SuperK. The widths of the shaded regions represent uncertainties in the average neutrino energy from 
the use of a diffusion scheme for neutrino transport (for details, see Refs. |275[I276] ). Observable effects 
of quarks only become apparent for stars older than 10 to 20 s. Sufficiently massive stars containing 
negatively charged, strongly interacting, particles (such as quarks, but also including hyperons and kaon 
condensates) may collapse to black holes during the first minute of evolution. Since the neutrino flux 
vanishes when a black hole forms, this would constitute an obvious signal that quarks (or other types 
of strange matter) have appeared. The collapse timescales for stars containing quarks are predicted to 
be intermediate between those containing hyperons and kaon condensates. 
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Figure 40: Total neutrino luminosity for proto neutron stars made of nucleons and quarks 
275 . Reprinted figure with permission from J. A. Pons et ai, Phys. Rev. Lett. 86 (2001) 
5223. Copyright 2001 by the American Physical Society. 

Already a few hours after birth, the internal neutron star temperature drops to ~ 10 9 K. The 
cooling of the star is primarily dependent for the next several thousand years on the neutrino emissivity 
of the core's composition. Photon emission overtakes neutrinos only when the internal temperature 
has fallen to ~ 10 8 K, with a corresponding surface temperature roughly two orders of magnitude 
smaller. Being sensitive to the adopted nuclear equation of state, the neutron star mass, the assumed 
magnetic field strength, the possible existence of superfluidity, meson condensates and quark matter, 
theoretical cooling calculations provide most valuable information about the interior matter and neutron 
star structure. The thermal evolution of a neutron star also yields information about such temperature 
sensitive properties as transport coefficients, transition to superfluid states, crust solidification, and 
internal pulsar heating mechanisms such as frictional dissipation at the crust-superfluid interfaces s 277j. 
In general, the possible existence of meson condensates and (certain phases of) quark matter enhance 
the neutrino emissivity from the core, leading to a more rapid early cooling. An overview of neutrino 
emitting processes from the core and crust of a neutron star is given in table |H1 Superfluidity of nucleons, 
on the other hand, has the opposite effect on cooling. Quantitative constraints on cooling have been 
hampered by the relatively small number of young pulsars known, the complication that several of 
them also display non-thermal, beamed X-ray emission from their magnetospheres, and uncertainties 
in distance and interstellar absorption. Table [7| summarizes the temperatures of a collection of neutron 
stars. To compare the observations with theory, one needs mainly the neutron star effective temperatures 
T s and ages r, which are compiled in table[7|for a representative collection of neutron stars. The thermal 
photon luminosity in the local reference frame of the star is given by L 7 = 4nR 2 aT^ with a the Stefan- 
Boltzmann constant. The apparent (redshifted) effective surface temperature T s °° and luminosity L^°, 
as detected by a distant observer, are given by 

T s °° = T s ^Jl- RJR and L~ = 47^°°) V(T & °°) 4 = L 7 (1 - R s /R), (116) 
where R s = 2M = 2.95 M/M Q km is the Schwarzschild radius. The coordinate radius R is connected 
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Table 6: Overview of neutrino emitting processes relevant for neutron star cooling [2|. 



Name 


Processes 


Emissivity 


Efficiency 


Modmed Urea 


n + n^n + p + e~~ + v e 
n + p + e~ ^n + n + ^ e 


~ 10 i g 


slow 


Direct Urea 


n — > p + e~ + P e 
p + e~ — > n + z/ e 


~ 10 27 T| 


fast 


Quark modified Urea 


a + u + e a + a + z/ e 
u + u + e~ —>■ u + d + u e 
d + u + e~ — > d + s + v e 
u + u + e~ — > m + s + v e 


~ 10 lg 


slow 


/ \ 1 1 • i TT 

Quark direct Urea 


1 i — i — 

d — > u + e + z/ e 

i — 7 i 

it + e — > a + z/ e 
s — * m + e~ + z/ e 
u + e~ — > s + z/ e 


~ 10 °T 9 


fast 


7r~ condensate 


n+ < 7r~ >^ n + e~ + z/ e 


~ 10 26 T 9 6 


fast 


K~ condensate 


n+ < ff~ >^ n + e~ + z/ e 


~ 10 26 T 9 6 


fast 


Quark bremsstrahlung 


Qi + V2 — ► Ql + Q2 + v + v 


~ 10 Jg 


slow 


Core bremsstrahlung 


n + n^n + n + u e + u e 
n + p — ► n + p + v e + v e 
e~ + p — > e~ + p + z/ e + v e 


~ 10 19 T 9 8 


slow 


Crust bremsstrahlung 


e~ + (A,Z) -+e- + (A,Z) 
+u e + v e 




slow 



to the apparent radius according to Eq. ()108|) . 

As already mentioned in section I5.2.4[ recent Chandra X-ray observations have identified pulsar 
PSR 0205+6449 at the center of the young Crab-like supernova remnant 3C58. Historical evidence 
suggests an association of the remnant with supernova SN 1181, which makes 3C58 younger than Crab 
(see tabled). The Chandra observation indicate that the thermal component must be very small, since 
the radiation is nearly completely fit by a power- law spectrum |215j . The temperature of a possible 
residual thermal component is thereby limited to an effective black-body value of T°° < 95 eV (surface 
temperature of < 6.0334) |215j which, as can be seen from Fig. falls well below predictions of 
standard cooling calculations |215j . As pointed out by Prakash et al. |274j . this upper limit can be fit 
with standard neutrino cooling (like n + n — > n + p + e~ + u e ) plus pair-breaking and formation, but 
the luminosity and ages of other neutron stars cannot be simultaneously fit using the same equation of 
state |298| 1299] which has the interesting consequence that more exotic, rapid cooling processes may exist 
in the core of PSR 0205+6449. Physical processes which would enable such a rapid drop in temperature 
range from the presence of meson condensates |300j, to quark matter [101 130H I302j . to the direct 
Urea process |303j . (For a very recent overview on neutron star cooling, see Ref. |298j .) 

Before the discovery of color superconductivity of quark matter, it was believed that depending on 
the density of electrons in quark matter, the temperature of quark stars could drop much more rapidly 
than for neutron stars 1304| 1305] . The density of electrons is crucial since the fast quark direct Urea 
processes (see table [HJ), 

d— >u + e~ + v e , u + e~ — >d + is e , s— > u + e~ + v e , u + e~ — > s + v e , (H7) 

are only possible if the electron Fermi momentum in quark matter is sufficiently high so that energy 
and momentum conservation in the above reactions is guaranteed. If the electron Fermi momentum is 
too small for this to happen, a bystander quark is needed to ensure energy and momentum conservation 
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Table 7: Surface temperatures of neutron stars at infinite distance from the stars, r denotes 
the age of the star 
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in the scattering process. The emissivity in the latter case is considerably smaller than the emissivities 
associated with the direct Urea processes in ()117j) . because of the different phase spaces associated with 
two-quark scattering and quark decay. If the electron fraction in quark matter vanishes entirely both 
the quark direct and the quark modified Urea processes become unimportant. The neutrino emission 
is then entirely dominated by bremsstrahlung, 

Q1 + Q2 — >Qi + Q 2 + v + v } (118) 

where Q\, Q2 denote any pair of quark flavors. In this case quark star cooling would proceed rather 
slowly, at about the same rate as conventional neutron stars ffi 1304} l305j . The same were the case if 
quark matter would be superfluid with gaps on the order of a few MeV, as described in the paper by 
Bailin and Love |306|I307] . The neutrino emissivities would then be suppress by a factor of exp(— A/T), 
with A the gap energy. 

The situation is more complicated if quark matter forms a color superconductor [2H1- If, as in the 
CFL phase, all quarks have a gap A ^> T then both the heat capacity Cy and neutrino luminosity L v 
are suppressed by ~ exp(— A/T) which would render quarks in the centers of compact stars invisible. 
Vanishingly small quark gaps, on the other hand, would lead to cooling behaviors indistinguishable 
from those of ordinary neutron stars made of either nucleons or nucleons and hyperons. In Ref. 308J, 
the photon and neutrino emission rates from the decay of photons and Nambu-Goldstone (NG) bosons 
associated with the spontaneous breaking of baryon number, U(1)b, in the CFL phase were calculated. 
The emission rates were found to be very small so that they would be inefficient for the core cooling 
of neutron stars containing quark matter in the CFL phase, rendering quark pairing in the CFL phase 
invisible to telescopes. This finding is in accordance with the quantitative determination of the mean free 
paths and thermal conductivities of photons (7) and NG bosons (0) in the CFL phase performed in 309J. 
The total conductivity associated with these particles was found to be kqfl = + k 7 ~ (27r 2 /9)T 3 _R 
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Figure 41: Cooling behavior of a 1.4 M & neutron star for competing assumptions about the 
properties of superdense matter. Three distinct cooling scenarios, referred to as 'standard', 
'intermediate', and 'enhanced' can be distinguished. The band-like structures reflects the 
uncertainties inherent in the equation of state of superdense matter [2|. 

which can be conveniently written as 

kcfl ^ 1-2 x 10 32 (T/MeV) 3 (R /km) erg cm" 1 sec" 1 K" 1 , (119) 

where Rq is the radius of the quark matter core. This expression reveals that the thermal conductivity 
of the CFL phase from photons and NG bosons is many orders of magnitude larger than the thermal 
conductivity of regular nuclear matter in a neutron star. The cooling of the quark matter core in the 
center of a compact star thus arises primarily from the heat flux across the surface of direct contact 
with the nuclear matter enveloping the CFL quark matter core. Since the thermal conductivity of the 
neighboring layer is also high, the entire interior of the star should be nearly isothermal |309| . The 
results in Ref. [309} confirm that the cooling time of neutron stars with CFL quark matter cores is 
similar to that of conventional neutron stars. 

Finally we mention briefly the cooling behavior of compact stars hiding color superconducting 2SC 
quark matter in their cores. The cooling of such stars is complicated by the fact that up and down quarks 
may pair with a gap A ~ 100 MeV that is orders of magnitude larger than the stellar temperature, 
£ 1 MeV, and are therefore inert with respect to the star's temperature evolution. In contrast to the 
CFL phase, where diquark condensation produces gaps for quarks of all three flavors and colors, there 
exits quark pairing channels that lead to weak pairing with gaps on the order of several keV to about 
1 MeV, which is on the same order of magnitude as the star's temperature. These quarks may thus not 
pair but, instead, may radiate neutrinos rapidly via the quark direct Urea process shown in table El If 
this is the case, the 2SC quark matter core will cool rapidly and determine the cooling history of the 
star [2H I302j . Examples of cooling curves of neutron stars containing quark matter in the 2SC phase 
are shown in Figs. 021 and 031 for different star masses. A quark gap of 1 MeV, as chosen in Fig. 021 
leads to too slow a cooling, while a reduced gap of 50 keV reproduces the observed data quite well. 
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Figure 42: Cooling curves of 2SC quark Figure 43: Same as Fig.l4*2lbut for a quark 

hybrid stars for a quark pairing gap of pairing gap of 50 keV. (From Ref. |3U2j .) 

1 MeV. (From Ref. [3712] .) 

7 Signals of Quark Matter in Rotating Neutron Stars 

In this section we explore possible signals of quark deconfmement in neutron stars, assuming that the 
densities in the centers of such objects are high enough so that quark deconfmement occurs. A convincing 
discovery of the type of signals described in this section could indicate that strange quark matter is not 
absolutely stable, ruling out the absolute stability of strange quark matter and the existence of strange 
quark stars, since it is impossible for (quark-hybrid) neutron stars and strange quark stars to both 
coexist stably. The signals described below require rather pronounced modifications of the equation 
of state caused by quark confinement and depend on the rate at which the mixed phase of quarks 
and hadrons gives way to pure quark matter. In addition great care is to be taken with respect to 
the numerical modeling |31()j of rotating stars as well as the properties of the nuclear curst |311j . We 
shall begin our discussion with isolated rotating neutron stars, which spin down because of the loss of 
rotational energy caused by the emission of an electron-positron wind from the star and by the emission 
of electromagnetic dipole radiation. This is followed by a discussion of accreting neutron stars in binary 
system. The spin period of such neutron stars increases over time. The densities inside both neutron 
stars that are spinning down as well as neutron stars that are being spun up by mass accretion changes 
dramatically, which could lead to signals of quark matter in observable data. 

7.1 Isolated Pulsars 

It is known from Figs. |^ and 021 that the weakening of centrifuge accompanied by the slowing down of 
a rotating neutron star causes a significant increase of its central density. From Fig. |^ for instance, 
one reads off that the central density of a neutron star model of mass M = 1.42 M Q , computed for a soft 
equation of state(Gg^gQ in the present case [312] ). increases from about 450 MeV/fm 3 for rotation at the 
mass shedding frequency, Ok, to more than 1500 MeV/fni for at zero rotation, which is a ~ 66% effect. 
Such dramatic changes in the interior density of a neutron star driven by changes in frequency modify 
the stellar composition considerably. If the mass and initial rotational frequency of a pulsar is such 
that during its slowing down phase the interior density rises from below to above the critical density 
for the quark-hadron phase transition, first at the center where the density is highest (Figs. HH1 EH1 
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Figure 44: Central star density versus ro- 
tational frequency for neutron stars of dif- 
ferent masses. The vertical bars mark the 
density where quark matter is produced. 
(From Ref. 0) 



n/n K 

Figure 45: Same as Fig. |EJ but for equa- 
tion of state Gg?g§ (From Ref. 0.) 



through l4T?j) and then in a region expanding in the radial outward direction away from the star's center, 
matter will be gradually converted from the relatively incompressible nuclear matter phase to the more 
compressible quark matter phase, as shown in Fig. |SU1 The tremendous weight of the overlaying layers 
of nuclear matter tends to compress the quark matter core, which causes the entire star to shrink on 
a length scale of several hundred meters, as shown in Figs. EB1 through H3 The mass concentration in 
the core will be further enhanced by the increasing gravitational attraction of the quark core on the 
overlaying nuclear matter. The moment of inertia thus decreases anomalously with decreasing rotational 
frequency as the new phase slowly engulfs a growing fraction of the star |312j . as can be seem from 
Figs. 03 and E21 Figure IB*T1 shows the moment of inertia, J, computed self-consistently from Eq. of 
several sample stars having the same baryon number but different internal constitutions [313J. The curve 
labeled M = 1.420 M©, computed for Ggf§Q, shows the moment of inertia of the quark-hybrid star of 
Figs. HH1 and The other curves correspond to a standard hyperon star (n, p, H) constructed for G^|° 
and a standard neutron star (n,p) where hyperons (H) have been ignored purposely. In accordance 
with what has been said just above, the shrinkage of quark-hybrid stars driven by the development of 
quark matter cores is the less pronounced the smaller the quark matter cores which are being built up 
in their centers. Correspondingly, the dip in I weakens with decreasing star mass, as shown in Fig. |^ 
for several sample masses in the range 1.416 < M/M Q < 1.420. Model calculations indicate that very 
strong reductions of J, such as found for the 1.421 M & model, for instance, may hardly be obtainable for 
physical scenarios other than the hypothetical quark-hadron phase transition 0. Hyperon populations 
alone as calculated in Ref. jT3] (Figs. EH ESI and E3J), for instance, appear to modify the equation 
of state way too little in order to cause significant changes in /. Nevertheless, there are models for 
the equation of state of hyperonic matter which also can strongly affect the spin evolution of isolated 
neutron stars |314t I315j . As shown in these references, depending on the nucleon-hyperon interaction 
and hyperon-hyperon interaction in matter, hyperons can even cause backbending. The backbending 
episode can terminate either unstably or through a stable continuous transition to a standard spin- 
down behavior. The observation of backbending in the timing behavior of isolated pulsars, or of spin 
clustering on accreting neutron stars discussed in section 17.21 is therefore not unambiguous evidence 
for quark deconfinement. It is also explored to which extent backbending may be caused by other 
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Figure 46: Frequency dependence of 
quark structures in equatorial star direc- 
tion for equation of state GgfgQ and a non- 
rotating star mass of 1.416 M & [2|. 
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direction [2] . 



competing particle processes (cf. Figs. IT1andl^T|). 

The decrease of the moment of inertia caused by the quark-hadron phase transition, shown in Fig. 1511 
is superimposed on the response of the stellar shape to a decreasing centrifugal force as the star spins 
down due to the loss of rotational energy. In order to conserve angular momentum not carried off 
by particle radiation from the star, the deceleration rate Cl (< 0) must respond correspondingly by 
decreasing in absolute magnitude. More than that, Cl may even change sign, as shown in Fig. E21 312J, 
which carries the important astrophysical information that an isolated pulsar may spin up during a 
certain period of its stellar evolution. The situation may be compared with an ice skater who spins up 
upon contraction of the arms before air resistance and friction of the skate on the ice reestablishes spin- 
down again. Such an anomalous decrease of / is analogous to the 'backbending' phenomenon known 
from nuclear physics, in which case the moment of inertia of an atomic nucleus changes anomalously 
because of a change in phase from a nucleon spin-aligned state at high angular momentum to a pair- 
correlated superfluid phase at low angular momentum. In the nuclear physics case, the backbending in 
the rotational bands of nuclei that was predicted by Mottelson and Valatin |316j and then observed years 
later by Stephens and Simon |317j . and Johnson, Ride and Hjorth [318J. For neutron stars, the stellar 
backbending of I is shown in Fig.|S21 Stars evolving from 'b' to 'a' are rotationally accelerated (Cl > 0), 
while stars evolving from 'a' to 'b', which could be part of the evolutionary track of pulsars accreting 
matter from companions (see section I7~2)) . are rotationally decelerated (Cl < 0). As we shall see next, 
the structure in the moment of inertia and, specifically, the backbending phenomenon dramatically 
modifies the timing structure of pulsar spin-down, rendering the observation of quark matter in neutron 
stars accessible to radio astronomy. Pulsars are identified by their periodic signal believed to be due to 
a strong magnetic field fixed in the star and oriented at an angle from the rotation axis. The period 
of the signal is therefore that of the rotation of the star. The angular velocity of rotation decreases 
slowly but measurably over time, and usually the first and occasionally second time derivative can also 
be measured. Various energy loss mechanisms could be at play such as the magnetic dipole radiation, 
part of which is detected on each revolution, as well as other losses such as ejection of charged particles 
319 . If one assume that pulsar slow-down is governed by a single mechanism, or several mechanisms 
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quark structure in equatorial neutron star 
direction computed for G^gg and a non- 
rotating star mass of 1.42 M [313J. 
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Figure 49: Same as Fig. HBJ but in polar 
direction |313j . 



having the same power law, the energy balance equation can then be written in the form 

f = s(5'< n > n, ) = - cnrH - (120) 

In the case of magnetic dipole radiation, the constant C is equal to C = |/i 2 sin 2 a where fi denotes the 
star's magnetic dipole moment. The quantity n in Eq. ()120|) is called braking index. It is n = 3 if / 
is kept constant during spin- up (down). If, as is customary, the star's angular velocity fl is regarded 
as the only time-dependent quantity, one obtains the usual formula for the rate of change of pulsar 
frequency, given by 

tt = -Ktt n , (121) 

with K = C 1 1 a constant. With the braking formula (j!21j) one can define the spin-down age of a pulsar 
given by 

T = -(n-iy^/Cl, (122) 

with n = 3 for energy loss governed by magnetic dipole radiation. However the moment of inertia is 
not constant in time but responds to changes in rotational frequency, as shown in Figs. [*""""| and E21 
more or less in accord with the softness or stiffness of the equation of state and according to whether 
the stellar mass is small or large. This response changes the value of the braking index in a frequency 
dependent manner, that is n = n(Q), even if the sole energy- loss mechanism were pure magnetic dipole, 
as expressed in Eq. (jl20j) . Thus during any epoch of observation, the braking index will be measured 
to be different from it canonical value n = 3 by a certain amount. How much less depends, for any 
given pulsar, on its rotational frequency and for different pulsars of the same frequency, on their mass 
and on their internal constitution 119U| I321[ I322[ I323j . When the frequency response of the moment 
of inertia is taken into account, Eq. (|121j) is replaced with [2*1 1312j 

ft = -2CQ n {21 + Qidl/dQ))' 1 . (123) 

This explicitly shows that the frequency dependence of Q corresponding to any mechanism that absorbs 
(or deposits) rotational energy cannot be a simple power law as given in Eq. (|121|) (with K a constant). 
It must depend on the mass and internal constitution of the star through the response of the moment 
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Figure 50: Fraction of mass existing in the form of pure quark matter, pure hadronic matter, 
and in the mixed phase of quarks and hadrons for the star of Fig. (From Ref. 0.) 

of inertia to rotation as expressed in ()123|) . Equation (|123|) can be represented in the form of (|121|) . but 
now with a frequency dependent prefactor, by evaluating 

n(n) - 2$ = 3 - ( 3 £ n + £L nA(2 1 + g n V 1 . (124) 

One sees that this braking index depends explicitly and implicitly on Q. This relation reduces to the 
canonical expression n = 3 only if I is independent of frequency, which may not the case, as seen above, 
if there are compositional changes driven by a varying star frequency. As an example, we show in Fig. 
ESI the variation of the braking index with frequency for two selected quark-hybrid stars of Figs. IBTl and 
1521 For illustrational purposes we assume dipole radiation. Because of the response of the moment of 
inertia to quark deconfinement, the braking index deviates dramatically from the canonical value n = 3 
at rotational frequencies where quark deconfinement leads to the built-up of pure quark matter cores 
in the centers of these stars. Such anomalies in n(Q) are not obtained for conventional neutron stars or 
hyperon stars because their moments of inertia increase smoothly with Q, as known from Fig. The 
observation of such an anomaly in the timing structure of pulsar spin-down could thus be interpreted 
as a signal of quark deconfinement in the centers of pulsars. Of course, because of the extremely small 
temporal change of a pulsar's rotational period, one cannot measure the shape of the curve which is 
in fact not necessary. Just a single anomalous value of n that differed significantly from the canonical 
value of n = 3 would suffice |3121 1324j . 

Carried over to the observed pulsar data of Q and Q, it appears that the change in centrifugal force 
over the life of a canonical, slowly-rotating pulsar could eventually be too meager to span a significant 
change. The significant braking anomaly, therefore, could be restricted to millisecond pulsars. For 
them, the phase change may occur only in such millisecond pulsars that rotate near the maximum-mass 
peak determined by the underlying equation of state. Otherwise the fraction of pure quark matter in 
their centers may not be sufficient to cause the required shrinkage. The phase change itself may be first 
(as in our example) or second order. Both orders will cause a signal as long as quark deconfinement 
causes a sufficient softening of the equation of state and quark matter is generated at the center of 
the star at a sufficiently high rate. On the observational side, a serious drawback may be that the 
braking indices of millisecond pulsars are very hard to measure, because of timing noise which renders 
the determination of Q very complicated. As a final but very important point on the subject of quark 
deconfinement, we estimate the typical duration over which the braking index is anomalous if quark 
deconfinement is well pronounced, as for the quark-hybrid star of mass M(0) = 1.421 M & . The time span 
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Figure 51: Moment of inertia versus ro- 
tational frequency of neutron stars having 
the same baryon number, A, but different 
internal constitutions. The dips at low fi's 
are caused by quark deconfinement [TQUj . 
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can be estimated from AT ~ — AQ/Cl = AP/P , where AQ is the frequency interval of the anomaly. 
The range over which n(O) is smaller than zero and larger than six (Fig. I53j) is Afl m — 100 s _1 , or 
AP ph —2irAQ/Q 2 rs 3 x 10~ 4 s at Q = 1370 s _1 . Hence, for a millisecond pulsar whose period derivative 
is typically P ~ 10 -19 one has AT ~ 10 8 years, as graphically illustrated in Fig. [SHI The dipole age of 
such pulsars is about 10 9 years. So, as a rough estimate we expect about 10% of the ~ 30 presently 
known solitary millisecond pulsars to be in the transition epoch during which pure quark matter cores 
are gradually being built up in their centers. These pulsars could be signaling the ongoing process of 
quark deconfinement in their cores. Last but not least we note that the spin-up time (region b-a in 
Fig. is about 1/5 of the time span AT, or about 1/50 of the dipole age. To avoid confusion, we 
point out that the spin-up has nothing to do with the minuscule spin-ups known as pulsar glitches. In 
the latter case the relative change of the moment of inertia is very small, AI/I ~ — AQ/Q ~ 10~ 6 
or smaller, and approximates closely a continuous response of the star to changing frequency on any 
time scale that is large compared to the glitch and recovery interval. Excursion of such a magnitude as 
quoted would fall within the thickness of the line in Fig. 

7.2 Accreting Neutron Stars 

The signal of quark deconfinement described in section 17.11 is computed for isolated neutron stars, 
where deconfinement is driven by the gradual stellar contraction as the star spins down. The situation 
is reversed in neutron stars in binary systems, which experience a spin-up torque because of the transfer 
of angular momentum carried by the matter picked up by the star's magnetic field from the surrounding 
accretion disk |320t I321| 1322} I325| I326j . The spin- up torque causes a change in the stars' angular 
momentum that can be expressed as |320j 

dJ/dt = Ml(r m ) - N(r c ) , (125) 

where M denotes the accretion rate and 

l(r m ) = ^[MrZ (126) 
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Figure 53: Frequency dependence of hy- 
peron thresholds in equatorial neutron 
star direction computed for HV. The star's 
non-rotating mass is 1.40 M Q [2|. 
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Figure 54: Same as Fig. E31 but in polar 
direction [2] . 



is the angular momentum added to the star per unit mass of accreted matter. The quantity N stands 
for the magnetic plus viscous torque term, 

N{r c ) = Kfi 2 r- 3 , (127) 

with [L = R 3 B the star's magnetic moment. The quantities r m and r c denote the radius of the inner 
edge of the accretion disk and the co-rotating radius, respectively, and are given by 

r m = £r A , r c =(Mfi- 2 ) 1/3 , (128) 
with (£ ~ 1). The Alfven radius r& is defined by 

r A = (fi 4 {2MM 2 y 1 ) 1/7 . (129) 

Accretion will be inhibited by a centrifugal barrier if the neutron star's magnetosphere rotates faster 
than the Kepler frequency at the magnetosphere. Hence r m < r c , otherwise accretion onto the star will 
cease. The rate of change of a star's angular frequency Q then follows from Eq. ()125j) as 



with the explicit time dependences as indicated. There are two terms on the right-hand-side of Eq. 
(jl30J) that grow linearly respectively quadratically with Q. Ignoring the linear term shows that mass 
transfer can spin up a neutron star to an equilibrium period of [230 

where Rq and B 9 are the star's radius and its magnetic field in units of 10 6 cm and 10 9 G, respectively. 
Medd in Eq. (J131J) denotes the maximum possible accretion rate, defined by the Eddington limit (see Eq. 
(|11U|) ). at which the accretion luminosity equals the luminosity at which the radiation pressure force on 



67 




Figure 55: Braking index, n, of quark- 
hybrid stars of Figs. ED and E21 The 
anomalies in n at Q ~ 400 s _1 and ~ 
1370 s" 1 are caused by quark deconfine- 
ment. The overall reduction of n below 3 
is due to rotation. (From Ref. [2].) 
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Figure 56: Braking index versus time 
for the quark-hybrid star of mass M = 
1.421 M & of Fig. ESJ The epoch over which 
n is anomalous because of quark decon- 
finement, ~ 10 8 years, is indicated by the 
shaded area. (From Ref. j2].) 



ionized hydrogen plasma near the star balances the gravitational acceleration force exerted by the star 
on the plasma. This condition leads to an Eddington accretion rate of M^dd = 1.5 x 10~ 8 i?6^© Y^ 1 - 
For a typical accretion rate of M_ 10 = M/(10 _10 M Q yr -1 ), the Eddington rate can be expressed as 
M E dd = 150 R 6 MZ\ M. The low -mass X-ray binaries (LMXBs) observed with the RXTE are divided 
into Z sources and A(toll) sources, which accrete at rates of M_io ~ 200 and M_ 10 ~ 2, respectively 

The solution of Eq. (|130|) in combination with the expression for the moment of inertia derived in 
Eq. (jHSl) for the quark-hybrid model Ggfgg (M(0) = 1.42 M ) is shown in Fig. E3 The magnetic field 
is assumed to evolve according to 

B(t) = B(oo) + (B(t = 0) - B(oo)) e~ t/td , (132) 

with t = at the start of accretion, B(t = 0) = 10 12 G, B(oo) = 10 s G, and t d = 10 6 yr. Such a 
decay to an asymptotic value seems to be a feature of some treatments of the magnetic field evolution 
of accreting neutron stars [327 . Moreover, it expresses the fact that canonical neutron stars have high 
magnetic fields, ~ 10 12 G, and millisecond pulsars have low fields of ~ 10 8 G. The result for the spin-up 
of the quark-hybrid stars is most striking. One sees that quark matter remains relatively dormant in 
the stellar core until the star has been spun up to frequencies at which the central density is about 
to drop below the threshold density at which quark matter is predicted to exist for this model. As 
known from Fig. E21 this manifests itself in a significant increase of the star's moment of inertia. The 
angular momentum added to a neutron star during this phase of evolution is therefore consumed by 
the star's expansion, inhibiting a further spin-up until the entire quark matter core has been spun out 
of the center, leaving the star with a mixed phase of quarks and hadrons made up of hadrons and 
quarks surrounded by ordinary nuclear matter (see Figs. I4*H1 and 149)) . Such accreters, therefore, tend 
to spend a much greater length of time in the critical frequencies than otherwise. There will be an 
anomalous number of accreters that appear near the same frequency, as shown in Fig. EH1 Evidence 
that accreting neutron stars pile up at certain frequencies, which are well below the mass shedding limit, 
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Figure 57: Evolution of spin frequencies of accreting X-ray neutron stars with (solid curves) 
and without (dashed curves) quark deconfinement. The spin plateau around 200 Hz signals 
the ongoing process of quark reconfinement in the stellar centers. (From Ref. |320j .) 

is provided by the spin distribution of accreting millisecond pulsars in 57 Tuc and neutron stars in low 
mass X-ray binaries observed with the Rossi X-ray Timing Explorer. The proposed limiting mechanisms 
responsible for this behavior could be gravity-wave emission caused by the r-mode instability, or a small 
stellar mass quadrupole moment t 328, 329| I330| . As shown here, quark reconfinement may be linked to 
this phenomenon as well jS2Dl E2H E23 CUE] 



8 Summary 

The tremendous pressures in the cores of neutron stars might be able to break neutrons, protons plus 
other hadronic constituents in the centers of neutron stars into their quark constituents, creating a 
new state of matter known as quark matter which is being sought at the most powerful colliders. If 
quark matter exists in the cores of neutron stars, it will be a color superconductor whose complex 
condensation pattern is likely to change with density inside the star. The exploration of the numerous 
astrophysical facets of (color superconducting) quark matter is therefore of uppermost importance and 
is pursued by physicists from different, complementing fields of physics. Their joint scientific efforts, 
surveyed in this review, provides most valuable information about the phase diagram of nuclear matter 
at high baryon number density but low temperature, which is not accessible to relativistic heavy ion 
collision experiments, and may ultimately provide us with a glimpse on the kind of matter that filled 
our universe just milliseconds after is was born. 
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